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ABSTRACT 

The Carnegie Supernova Project (CSP) is designed to measure the luminosity distance for Type la 
supernovae (SNe la) as a function of redshift, and to set observational constraints on the dark energy 
contribution to the total energy content of the Universe. The CSP differs from other projects to date in 
its goal of providing an /-band rest-frame Hubble diagram. Here we present the first results from near- 
infrared (NIR) observations obtained using the Magellan Baadc telescope for SNe la with 0.1 < z < 0.7. 
We combine these results with those from the low-redshift CSP at z < 0.1 (Folatclli ct al. 200j|). In this 
paper, we describe the overall goals of this long-term program, the observing strategy, data reduction 
procedures, and treatment of systematic uncertainties. We present light curves and an /-band Hubble 
diagram for this first sample of 35 SNe la and we compare these data to 21 new SNe la at low redshift. 
These data support the conclusion that the expansion of the Universe is accelerating . When combined 
with independent results from baryon acoustic oscillations (jEisenstein et al.ll2005| ). these data yield 
0,„ = 0.27±0.02 (statistical), and Ude = 0.76±0.13 (statistical) ±0.09 (systematic), for the matter 
and dark energy densities, respectively. If we parameterize the data in terms of an equation of state, 
w (with no time dependence), assume a flat geometry, and combine with baryon acoustic oscillations, 
we find that w = —1.05 ± 0.13 (statistical) ± 0.09 (systematic). The largest source of systematic 
uncertainty on w arises from uncertainties in the photometric calibration, signaling the importance of 
securing more accurate photometric calibrations for future supernova cosmology programs. Finally, 
we conclude that either the dust affecting the luminosities of SNe la has a different extinction law 
(Rv = 1.8) than that in the Milky Way (where Ry = 3.1), or that there is an additional intrinsic color 
term with luminosity for SNe la, independent of the decline rate. Understanding and disentangling 
these effects is critical for minimizing the systematic uncertainties in future SN la cosmology studies. 
Subject headings: cosmology: observations - cosmology: distance scale - supernovae: general 
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1. INTRODUCTION 

Observations of high-rcdshift Type la super- 
novae (SNe la) currently provide the best e vi dence 
for an accelerating universe (|Riess et al.1 119981 : 
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Perlmutter et all Il999t iKnop et all 120031 ; lAstier et all 
20061 IWood-Vasev et all 120071 ). Independently, the 
Wilkinson Microwave Anisotropy Probe (WMAP) 
measurements of fluct uations in the cosm i c micr owave 
background (CMB) ((Spergel et all 120031 12001 : de- 
tections of acoustic oscillations in the matter density 
spectrum (E isenstein et al.l I2005D ; as well as the com- 
parison of Hubble expansion a ges from the Hubble Key 
Project (|Freedman et al. | I2001D with the ages of Milky 
Way globular clusters ( Krauss fc Chaboverll2003f ). have 



all led to the growing body of evidence for a dark-energy 
component to the overall mass-e n ergy density of the 
Universe (see also iFrieman et al.l (|2008l l for a recent 
review). The above studies yield results consistent with 
a flat universe where the sum of the dark energy and 
matter densities fljjE + = lj an< i where £Ide ~ 0.7, 
with a matter density f2 m ^0.3. 

At present, a physical understanding of this dark en- 
ergy component remains elusive, and there is a wide va- 
riety of possible alternatives ranging from the cosmolog- 
ical constant, as originally proposed by Einstein in 1917; 
a decaying scalar field; or perhaps even a modification of 
general relativity. A convenient measure is w, the ratio 
of pressure P to energy density p, where w = P j p. In 
the case of a cosmological constant, w = — 1. The time 
evolution of w is not known at present. Given our cur- 
rently limited understanding of this critical component 
of the Universe, it is clear that further observational and 
experimental data are needed to constrain and charac- 
terize the properties of dark energy. 

SNe la show a relation between peak brightness and 
rate of decline, in addition to color (or reddening). Ac- 
counting for these correlations yields a dispe rsion in 
the Hubble diagram of about 7% in distance ([Phillips 
19931 lHamuv et al.l|1995tlRiess et ai]|1996t lHamuv et al. 
1996a! lAstier et al.ll2006l ). As more and more SNe la 
are discovered and the statistical uncertainties in this 
method arc decreased, the challenge becomes under- 
standing and controlling the remaining systematic un- 
certainties, as well as testing for currently unknown sys- 
tematic effects. A well-known systematic is reddening 
and extinction due to dust, and there may potentially 
be differences due to chemical composition and age of 
the progenitors or their environment. The requirement 
for increasing measurement accuracy, the lack of a de- 
tailed theoretical understanding of SN la, the fact that 
most observations have been made at rest-frame opti- 
cal and ultraviolet colors (where reddening uncertainties 
are large), the difficulty of obtaining accurate K- and 
spectral-corrections, all point to the need to characterize 
and decrease these systematic errors and uncertainties. 
As other errors have been reduced over time, what were 
relatively small effects previously have now become in- 
creasingly important as the required precision in cosmol- 
ogy has risen. 

Observations of SNe la at near-infrared (NIR) wave- 
lengths offer several advantages in minimizing a number 
of systematic effects (most notably reddening), as well 
as providing an additional, independent probe of the ex- 
pansion history. We make use of the Y band, centered 
near 1.035 (im (|Hillenbrand et all I2002D . which falls in 
between the traditionally classified optical and infrared 
spectral regimes. The Y and J (1.25 fim) bands pro- 
vide a means of obtaining rest-frame /-band magnitudes 



for SNe la in the redshift range 0.1 to 0.7, an interval 
covering the time at which the influence of dark energy 
begins to dominate the expansion. Hence, NIR observa- 
tions offer an important opportunity both to minimize 
systematic uncertainties and yield an independent Hub- 
ble diagram. 

Aside from a few o bservations publish ed by 
iRiess et al. I (|2000D and iNobili et all (|2005f ). rest- 
frame /-band measurements have not yet been routinely 
undertaken at higher rcdshifts because at z^0.25 the 
/-band is rcdshiftcd beyond the CCD sensitivity limit at 
red wavelengths. A further challenge to NIR photometry 
obtained from the ground is that the sky background 
is much greater than at optical wavelengths. As a 
consequence, observing distant SNe la in the NIR is 
expensive in telescope time. However, since current 
searches to find SNe la now yield extensive optical 
coverage of the light curves, the decline rates and time 
of maximum light are already well defined, so that 
fewer observations are required to derive NIR light 
curve parameters accurately. With the availability of 
6.5-meter class telescopes and NIR arrays, a NIR study 
of SNe la at z > 0.1 has now become feasible. 

Ongoing optical surveys aimed at discovering large 
samples of SNe la are signficantly decreasing the sta- 
tistical uncertainties in SN la cosmology measurements, 
while simultaneously aiming to reduce the systematic er- 
rors. The CSP is complementary to these studies; that 
is, it is not designed to rival optical surveys on numbers 
of objects and statistical errors, but rather it is more 
focused on follow-up observations critical to testing for 
and minimizing systematic uncertainties. These follow- 
up observations (multiple wavelengths at low redshift, 
and NIR observations at higher redshifts) require more 
observing time per individual SN la. An understanding 
of the systematic errors is a critical goal for all SN la 
programs. 

To date, as part of the high-redshift CSP, we have ob- 
tained NIR observations of 80 SNe la. Subsequent host- 
galaxy observations have been obtained for 71 objects 
and we have fully analyzed 35 objects, the sample for 
which photometry, light curves and a preliminary Hub- 
ble diagram are presented in this paper. 

2. THE CARNEGIE SUPERNOVA PROJECT 

The CSP is an umbrella name for our two coordinated 
SN programs: one being carried out at the Las Cam- 
panas 1 m Swope and 2.5 m Dupont telescopes for low 
(0<z<0.1) redshift objects, and the other using the Mag- 
ellan Baade 6.5 m telescope at higher (0.1<z<0.7) red- 
shifts. Preliminary results from th e low-redshift part o f 
this project have bee n repo rted in lHamuv et alJ I2006h . 
IFolatelli et all (|2006l ) and iPhillips et al. I (|2007f) : the 
most recent results are presented in two compan ion pa- 
pers (|Contreras et al.l 12009b IFolatelli et all 12009( 1 . The 
CSP is not a SN search program: rath er as described be- 
low, as well as in lHamuv et al.l (|2006l ). we are using the 
facilities of the Las Campanas Observatory, in coordina- 
tion with several on-going search programs elsewhere, to 
provide detailed and immediate follow-up observations of 
newly-discovered SNe. 

2.1. Goals 
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Overall, the CSP is focused on obtaining observations 
of SNe la falling in the r edshif t range 0<z<0.7. As de- 
scribed in iHamuv et all (|2006[ ), at low redshifts we are 
currently obtaining 10-filter (uBV griY J H K s ) photome- 
try with excellent time sampling and optical spectroscopy 
to provide a database for the determination of the Hubble 
constant, accurate K-corrcctions, comparison with theo- 
retical models of SNe la, and a fundamental dataset for 
comparison with high redshift. At high redshifts, as de- 
scribed in this current paper, we arc obtaining Y J data 
near maximum light. The NIR data, in conjunction with 
optical photometry (being obtained as part of the search 
programs), allow the determination of reddening correc- 
tions and a rest-frame /-band Hubble diagram. One of 
the key goals of the CSP is to minimize the effects of 
reddening in the Hubble diagram, and ensure that the 
rest-frame (BVi) bandpasses, being observed at low red- 
shift, match those for a sample at higher redshift, so 
that reddening corrections can be applied in a uniform 
way. Ultimately the goal is to measure accurate lumi- 
nosity distances to better constrain cosmological mod- 
els. Assuming a model including dark energy, the goal 
is to characterize the acceleration of the Universe and 
constrain the equation of state, u>, to a precision and 
accuracy of ±10%. 

2.2. "I-band" filter nomenclature 

There are many "I-band" filters in current use. We 
distinguish here between the Kron-Cousins /-band, the 
SNLS McgaCam «M-band (based on Landolt standards 
calibrated to Vega), and the SDSS- II i<;-band (based o n 
a calibration tied to four subdwarfs. lSmith et al.l (|2002f )'). 
Our Las Campanas natural i-band calibration is also tied 
to Smith et al. The central wavelengths of these (and 
other) passbands are defined in Table [T] In this paper, 
we also refer to a generic 'I-band' when speaking collec- 
tively about observations in the NIR wavelength range 
~ 8000A. 

2.3. I-band Light Curves SNe la 

The form of NIR light curves for SNe la differs from 
those at optical wavelengths. At IYJH and K, the 
light curves show a minimum following maximum light, 
followed by an additional secondary maximum, less lu- 
minous than the primary peak (e.g.. lElias et al.l | 1985 |: 
Hamuv et al]|1996bl ; lMeikldl2000l ; iKrisciunas et allbooil ; 
Nobili et al.ll2005f l. This secondary maximum is thought 
to be the resul t of a sudden drop in th e mean opaci ty in 
the SN ejecta (jPinto fc Eastman! 120001 ; iKasenl 1200(1 . In 
Figure [1] we show some examples of i-band light curves 
for four well-observed, nearby SNe la from the low-z CSP 
sample. Our strategy for the high-redshift CSP is to ob- 
tain i photometry covering the first peak only, and to 
avoid the (fainter) secondary maximum. The optical sur- 
veys already provide measurements of the decline rate in 
the B band, so duplication of these data is not required. 
The median start and end times of our observations are 
and +12 days with respect to observed i-band maxi- 
mum, respectively. In this paper, unless otherwise noted, 
we refer to the observed (not rest) frame time. 

2.4. Targets 

There have been four on-going SN search programs 
with which the CSP has been actively collaborating. At 



low redshifts (z < 0.1), we are collaborating with the Lick 
Observatory S upern ova Search (LOSS), as described in 
IHamuv et alj (|2006f ). For intermediate redshifts (0.1 < z 
< 0.4), we have been collaborating with the Sloan Digital 
Sky Survey II (SDSSII), while the CFHT Legacy Survey 
(SNLS) and ESSENCE found SNe la out to higher red- 
shifts (0.1 < z < 0.7). The three surveys relevant to this 
highcr-rcdshift study are described briefly below. 

2.4.1. Sloan Digital Sky Survey (SDSS) II 

An extension t o the Sloan Dig ital Sky Survey 
(lYork et al. I 12 000). the SDSS-II, 23 (jSako et all [20081 : 
iFrieman et al.ll2~008h has completed a three- year rolling 
search program using the 2.5 m SDSS telescope to dis- 
cover intermediate redshift (0.05 < z < 0.35) SNe la in 
a 300 square degree southern equatorial field, obtaining 
multicolor (ugriz) photometry with a cadence of ~2- 
5 days. The data were obtained during the months of 
September through November 2005-2007. Spectra to de- 
termine redshifts and types for the SDSSII candidates 
were obtained with several telescopes with a range of 
apertures. Photometry for the 130 spectroscopically con- 
fir med SNe la from t he 20 05 season have been presented 
in lHoltzman et al. I (|2008| ). 

2.4.2. CFHT Legacy Survey (SNLS) 

The Supernova Legacy Survey 24 , a Canadian/French 
collaboration, used the CFHT and the 1 square-degree 
camera (MegaCam) to obtain deep optical (u'g'r'i'z') im- 
ages for four fields, each o f approximately on e square de- 
gree around the equator (jAstier et al.l l2006h . beginning 
in February, 2003. The SNLS was also a rolling search 
program in which each field was revisited every second 
night during a 5-month campaign each semester for a to- 
tal of five years. Spectroscopic follow-up for the SNLS 
candidates was obtained at the 10 m Keck telescope, the 
8 m Very Large Telescope (VLT) , the 8 m Gemini North 
and South telescopes, as well as the Magellan telescopes. 
In their first year of operation, they discovered 91 SNe la, 
and published data for a sample of 71 SNe la with red- 
shifts 0.1 < z < 1.1. In total, about 450 SNe la have 
been discovered. 

2.4.3. ESSENCE 

ESSENCE 25 completed its survey using the 4 m CTIO 
Blanco telescope and MOSAIC II wide-field camera cov- 
ering equatorial fields at VRI waveleng ths over the red- 
shift range between 15 < z < 0.75 (jMatheson et al.1 
120051 iMiknaitis et al.ll2007l ). The goal was to provide 
optical light curves for about 200 Type la SNe la over 5 
years (2002-2007). ESSENCE was scheduled for 30 half 
nights during a 3-month (October - December) campaign 
each year. The observing was centered on new moon, 
and was scheduled every other night, for a span of 20 
days, with a gap of 10 bright nights per month. Spec- 
troscopic follow-up for the ESSENCE candidates was 
obtained at Gemini North and South, Keck, the VLT, 
MMT, and Magellan telescopes. At the time of writ- 
ing, ESSENCE has publis hed observations for 102 SNe la 
(|Wood-Vasev et al.ll200l . 

23 http://www.sdss.org 

24 http:/ /www. cfht.hawaii.cdu/Scicncc/CFHLS 

25 http:/ /www. ctio.noao.edu/essence 
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2.5. CSP Magellan Follow-up Observations 

Our follow-up Magellan observations consist of Y- 
and/or J-band images. We chose not to observe at 
-ff-band because of the increased sky background pro- 
duced by atmospheric OH airglow emission. Optical pho- 
tometry for this sample of SNe la was obtained during 
the course of the original surveys (SDSS-II, ESSENCE, 
SNLS), and was not repeated. The Sloan i filter, which 
we are using at the Swope telescope for the low redshift 
part of the CSP, overlaps well with the Y-band for a 
redshift of z~0.35, and the J-band at redshift z~0.63 
(see Figure [5]) . If targets at these redshifts were avail- 
able, they moved to the top of our observing priority list. 
We aimed to obtain Magellan observations no more than 
a few days after maximum light, and preferably while 
still on the rise. Our criteria for selection of the SNLS, 
ESSENCE and SDSS-II objects were: 1) the SNe la were 
in the redshift range 0.2 < z < 0.7, 2) the SN la was es- 
timated to be near maximum, 3) SNe la near z = 0.35 
and z = 0.65 were given preference and 4) preference was 
given to SNe la well-separated from the host. In practice, 
criteria 3) and 4) were rarely invoked since there was not 
much choice after criteria 1) and 2). Our selection cri- 
teria therefore matched fairly closely those of the parent 
surveys. 

The SNLS, ESSENCE and SDSS-II projects provided 
coordinates, finder charts and epochs of high-probability 
SN la candidates for the Magellan observations. In three 
cases (SNLS 03Dlgl, SNLS 04D2cw, ESSENCE 2004-H- 
280), the spectral classifications were later revised (to 
types other than SNe la). In one case (SNLS 04D2al), 
the redshift was too large (z = 0.838) and the SN la was 
undetectable in our observations. In Table [H we list the 
SN la name, position, number of CSP observations, and 
IAU cross-identification, for the SNe la. 

3. INSTRUMENTATION AND OBSERVATIONS 

3.1. Instrumentation 

Our NIR imaging has been carried out using Y- and J- 
band filters on the Persson Auxiliary N asmyth Infrared 
Camera (PANIC) ([Martini et al.ll2004D . This camera is 
mounted on the east Nasmyth platform of the Magel- 
lan Baade telescope. It contains a 1024x1024 HgCdTe 
Hawaii-1 array with a scale of 0.125 arcsec pixel -1 , and 
covers a field of view of 2.1 'x 2.1' on the sky. 

3.2. Observations 
3.2.1. Observing Strategy and Current Status 

Our observing goal was to acquire Y and/or J photom- 
etry with gaps in the SN la light curves no larger than 
7 days, straddling the light curve over maximum light if 
possible. To avoid large gaps, PANIC was scheduled dur- 
ing dark/grey time, in addition to the customary bright 
time for infrared instruments. In practice, telescope time 
was scheduled with gaps generally less than 5 days, al- 
though gaps of 8 days were occasionally unavoidable. A 
record of observations, including finding charts of the ob- 
jects, is maintain ed on our web site for the pro ject, which 
can be found at |http: / / www.ociw.edu/cs p / . Generally, 
3 to 5 observations per SIN la were obtained. 

We began this long-term project with a pilot program 
in November 2003 running through April 2004, centered 
on the best weather period at Las Campanas. Our first 



host-galaxy images (once the SN la had faded) were ob- 
tained in the following year. In the second year, we ob- 
served from October 2004 through March 2005. Unfortu- 
nately, poor weather at Mauna Kea during the northern 
winters of 2003 and 2004 resulted in a very low yield of 
SNe la for CSP follow-up during this pilot project. More- 
over, ESSENCE was operational only during the north- 
ern fall, and SDSS-II began routine operations during 
the fall of 2005. We therefore shifted the CSP follow-up 
campaigns to August through January when there was 
good overlap with all three SN search programs. 

After four campaigns, we obtained Y and/or J photom- 
etry for a total of 80 objects, 75 of which were ultimately 
confirmed as SNe la. Follow-up host-galaxy observations 
have been obtained for 71 of these SNe la and we are in 
the process of obtaining host-galaxy images for the re- 
maining objects. The redshift distribution for our total 
sample is shown in Figure [3j In this paper, we report on 
the photometry from the first 35 SNe la. A summary of 
the observations for these 35 SNe la is given in Table [51 
At the current time, data for 20, 13, and 2 SNe la have 
been analyzed from the SNLS, SDSS-II and ESSENCE 
surveys, respectively The original goal of the CSP was 
to acquire a sample of 100 objects over the redshift range 
0.1 < z < 0.7; however, with the conclusion of the three 
surveys, our expected sample will be limited to 75 SNe la. 

3.2.2. Observing Procedures 

Twilight sky flats were taken in each filter every night 
in at least 5 dithered positions with 2 exposures at each 
position, controlling the count level to be in the linear 
regime. The detector becomes nonlinear at the 1% level 
at approximately 13,000 ADU and at the 5% level at 
approximately 35,000 ADU. Dark frames were taken in 
sets of 15, each with exposure times matching those of the 
science images. The SN la observations were obtained 
using 9 dither positions with 2 exposures per position. 
The exposure times used for individual images were 120 
seconds with total exposure times ranging from 2160 to 
8640 sec. Reference images for these SN la fields were 
taken in the same manner once the SN la faded, to allow 
accurate background subtraction. The signal to noise 
in the stacked reference images were as high (or higher) 
than those in the SN la images, so as not to degrade the 
photometry. The seeing for these images ranged typically 
between 0.4 and 0.8 arcsec. On photometric nights, 3 
to 5 stars chosen fr om the standard system defined by 
iPersson et al.l (|1998h were observed throughout the night 
in each filter. These standard stars were observed in 5 
dither positions with 2 exposures at each position. The 
exposure times of individual images were 3 or 4 seconds, 
resulting in total exposure times of 30 and 40 seconds. 
On occasion, the telescope was de-focused to ensure the 
counts remained in the linear regime of the detector. 

4. DATA ANALYSIS 

4.1. Pipeline Processing 

We have developed an automated pipeline to process 
the PANIC images. The pipeline carries out (1) linear- 
ity corrections, (2) dark combination and subtraction, 
(3) bad-pixel mask production, (4) flat-field combination 
and division, (5) sky image computation and subtrac- 
tion and (6) combination of dithered frames into final 
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stacked images. We apply a predetermined linearity cor- 
rection law to every pixel value above 8,000 ADU. The 
multiplicative correction ranges from 1.0 at 8000 ADU 
to 1.06 at 40,000 ADU. Dome flats are created by sub- 
tracting images of equal exposure times taken with no 
dome lamps from those with the dome lamps on. The fi- 
nal dome flats are used only to build bad pixel masks for 
each night and median-combined twilight sky flats are 
used for flat-fielding. Sky frames are subtracted from 
the individual object frames using modal scale factors. 
Finally, stacked images are created by aligning and aver- 
aging the individual object frames. 

4.2. Galaxy Template Subtraction 

To obtain accurate photometry for the SNe la, the host 
galaxy light must be subtracted from the images. The 
strategy for removing the host galaxies from the SN la 
data involves three steps: (1) obtaining host-galaxy im- 
ages of the SN la fields in the year following the events; 
(2) registering and matching the point-spread functions 
(PSF) of the images; (3) subtracting the PSF-matched 
SN la and host-galaxy images. The algorithm developed 
for the registration and non-parametric matching of the 
PSFs will be descr ibed in more detail in a later paper 
([Kelson et alJ 12009); the method is summarized briefly 
here. 

For a given SN la image and associated host-galaxy im- 
age, SExtractor (jBertin fc Arnoutslll996h is used to iden- 
tify objects down to a threshold of 3-er. The positions of 
these objects are used to compute the coordinate trans- 
formation between each SN la image and the host-galaxy 
image. Matching the PSF of the host-galaxy image, T, 
to that of the SN la image, S, consists of solving for the 
convolution kernel, k, that maps point-sources in T to 
S. Once k is known, the entire host-galaxy image can 
then be convolved with k and subtracted from S to iso- 
late the SN la. However, our fields of view are small and 
in several instances lack suitable point sources for de- 
termining the kernel using more traditional techniques, 
for example, parameterizing the kernel with a Gaussian 
(Alard & Lupton 1998, Alard 2000). Instead, we use a 
non-parametric technique that utilizes all objects in the 
field to constrain the kernel. The kernel consists of a 
(2M + 1) x (2M + 1) matrix indexed by u and v, and is 
determined by minimizing 



N 

E 



S{xi,yi)- £„ k (u, v) [T(a 



v)] 



a{xi,Vi) 



where the sum is over N pixels in the source and template 
images. 

The matrix of T(xi~ u, Di—v) is decomposed using sin- 
gular value decomposition, and its constituent eigenvec- 
tors contain the kernel's natural set of orthogonal basis 
functions. The eigenvalues represent each eigenvector's 
sensitivity to noise. We eliminate those basis functions 
that do not contribute to reducing the \ 2 P er degree of 
freedom. Because image re-binning is mathematically 
equivalent to a convolution, errors in the registration are 
fully accounted for by the convolution kernel that mini- 
mizes x 2 - 

Ideally, one wishes to convolve (i.e., degrade) the host- 
galaxy images so as not to decrease the signal-to-noise of 



the SN la. However, approximately 8% of our SN la im- 
ages were taken under exceptional seeing conditions (less 
than 0.35 arcsec) which, to date, have not been matched 
in our host-galaxy observations. As a result, this sub- 
set the SN la images have to be degraded to match the 
poorer image quality of the host-galaxy images. Fortu- 
nately, these SNe la have lower rcdshifts [z < 0.3) and 
are relatively bright, so that the signal to noise remains 
high. The SN la images in our sample to date have a 
range of image quality, from 0.3 arcsec to 1.4 arcsec; the 
host-galaxy images have a range of image quality, from 
0.3 arcsec to 1.3 arcsec. We find that our image subtrac- 
tion technique works well for removing the host galax- 
ies for a wide range of seeing conditions, position of the 
SN la relative to the host galaxy, and redshift. Three ex- 
amples of a host-galaxy-subtracted image are shown in 
Figure H for objects at redshifts of 0.25, 0.30, and 0.68, 
respectively. 

4.3. Photometry 

Observing SNe la in the NIR is more challenging than 
in the optical owing to higher sky background. The con- 
trast of SN la to galaxy is also less in the red, and over- 
coming these effects require longer integrations than in 
the optical. We measure the flux with a two-step ap- 
proach. The first step is to measure the magnitudes 
of several stars in each SN la field, which we shall re- 
fer to as tertiary standards. The pri mary standards ar e 
those that establish the JHK syste m (|Elias et ahl 1982D . 
on which the secondary standards (|Persson et al.Tl998f) 
that we observe are based. In choosing the tertiary stan- 
dards, we require that: 1) the star is no closer than 20 
arcsec to the edge of the PANIC field of view 2) the star 
appears in all observations of the SN la and the host 
galaxy; and 3) there are no significant residuals for the 
star after host-gal axy subtractio n. We use DAOPHOT 
and D AO GROW (jStetsonl [1990T ) to measure and fit a 
growth curve (flux versus aperture size) for each star. 
We then compute an aperture correction using these fits 
to get the flux measured through a 10-arcscc-diamctcr 
aperture, matching the aperture used on the standards. 
These fluxes are therefore calibrated with respect to the 
secondary standards observed that evening, and are then 
averaged over all photometric nights. 

The second step is to measure the flux ratios between 
the SN la and the tertiary standards. At the highest red- 
shifts, we are working at the detection limit of the tele- 
scope where DAOPHOT and D AO GROW are no longer 
robust. Instead, we u se the optimized extraction algo- 
rithm of (|Navlorlll996f ). In brief, the brighter stars in the 
field are used to estimate the PSF, P^j, modeled as a 
superposition of Gaussian and Moffat profiles. This es- 
timated PSF is then used as a weight mask in summing 
the flux from the SN la pixels: F = J2i w i,j ~ 



where 



P is the model PSF, D is the mea- 



sured counts, and S is the sky flux. Following iNavlorl 
(fl996h . we have assumed that the variance in the weights 
is dominated by the sky (Vij ~ V s in his equation (10)). 
This yields an optimized estimate of the flux for the 
SN la. We then use the same weight mask to measure 
fluxes and compute flux ratios between the SN la and 
each tertiary standard. These flux ratios, together with 
the calibrated fluxes from step 1, provide estimates of 
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the calibrated flux of the SN la. These estimates are av- 
eraged to yield the final calibrated flux for the SN la. We 
have done a number of tests in which fake point sources 
are inserted into the frames and then recovered, and find 
that the meth od works extremely w ell. These results will 
be reported in lKelson et alJ ()2009l ). 

We have also taken great care in computing the vari- 
ances of our fluxes. In the process of correcting for image 
distortion in the PANIC pipeline, rectifying the images, 
convolving with a kernel, and then subtracting a host 
galaxy image, the pixels in each image have become cor- 
related. We produce variance maps for each observed 
field early in the PANIC pipeline. A simple propagation 
of errors is done at each step of the pipeline, updating 
the variance in each pixel in the maps. The final vari- 
ance maps are used to compute the noise rather than 
image statistics, ensuring that we can properly estimate 
the variance of the tertiary standards and the SN la. 

4.4. Absolute Calibration 

The absolute calibration of the M agellan photometry i s 
based on JHK standard stars from lPersson et all (|1998l ) , 
where the zero point is tied to Vega. Currently, we are 
applying the K-ban d absolute calibration described in 
Hamuv et all (120061). using Kurucz model spectra and 
Hillenbrand et all (120021) standard star measurements. 



Contreras et al. ni2009Tl have recently obtained new ob- 



servations confirming this calibration to an accuracy of 
±0.01 mag. We have observed common standard stars 
and adopted the same procedures for the reduction of 
the standard stars so that differences in the calibration 
of low- and high-z samples are minimized. Adopting at- 
mospheric extinction coefficients of ky — 0.10 and kj = 
0.12 mag/airmass following Hamuy et al., we solve for 
the nightly zero points in each filter. Photometry for the 
standard stars is obtained with apertures of diameter 10 
arcsec. The statistical errors in the zero point are deter- 
mined from the scatter of individual measurements and 
range from ± 0.01 to ± 0.05 mag. A minimum of three 
photometric nights determines the absolute flux for each 
SN la (see M4.3[) . The average scatter in the zero points 
based on observations of standard stars is 0.023 mag at 
J and 0.013 mag at Y. 

In this (and other CSP) papers, we are presenting our 
photometry in our own natural system. The advantage 
of using the natural system is that it avoids the un- 
certainties resulting from the broad features present in 
the spectra of SNe la, which present challenges for the 
tr ansformation onto the standard system. As outlined 
in lContreras et all ()2009f ). the natural photometry is ob- 
tained by first computing color terms that transform lo- 
cal standard seque nces of sta r s to t he system appropri- 
ate for each filter (|Landolt I (|1992f ) stan dards for BV 
ISmith et alJ (|2002f ) standards for ugri, and lPersson et al.l 
(|1998t ) for YJ). These color terms are then used in re- 
verse to transform the standard magnitudes to our nat- 
ural system and it is these magnitudes that are used 
to calibrate the SN la photometry. The natural mag- 
nitudes are therefore equal to the standard magnitudes 
at zero color. Our natural system magnitudes can be 
straightforwardly transformed to other systems. Trans- 
mission curves for o ur N IR and optical fil t ers ar e given in 
lHamuv et all (|2006f ) and lContreras et all (|2009h . and up- 
dated versions are available online at the CSP web site. 



In the SN la target fields, we identified several isolated 
stars to serve as tertiary standards to determine absolute 
flux for those nights that were not photometric. Photom- 
etry for these secondary standards is also available on 
line. In Tabled! we list the filters, effective wavelengths, 
and published references to the calibrations relevant for 
the CSP. A comparison of our own NIR YJ photometry 
(K-corrected to the i band) to that of the optical surveys 
(SNLS, ESSENCE and SDSS-II) requires that we adopt 
a standard with a measured spectral energy distribution 
(SED) spanning th e optical and NIR. For th is purpose, 
we have chosen the Bohlin & Gilliland (2004) model for 
Vega, as updated in lBohlinl (|2007fh 



5. SYSTEMATIC EFFECTS 

As the number of objects has increased, SN la cosmol- 
ogy has reached the stage where the systematic uncer- 
tainties are becoming the dominant source of error. The 
CSP has been designed to minimize known systematic ef- 
fects, particularly those due to K-corrections, reddening, 
and cross calibrations to different photometric systems. 
Below we discuss our current approach to dealing with 
K-corrections and reddening. 

5.1. K-corrections 

The observed SEDs of SNe la are shifted and stretched 
with redshift due to the expansion of the Universe. Ac- 
curate corrections for these effects (K-corrections) for 
SNe la remain an observational challenge, and much ef- 
fort has been put into creating libraries of SN la spectra 
with which to assemble SEDs that can be used to esti- 
mate the K-corrections a nd their variance (| Nugent et al.l 
I2002t iHsiao et all 120071) . The K-correctio n s are com- 
puted in a manne r similar to iKim et all (jl996t ) and 
iNugent et ail (|2002|) . using the following formula: 



K AB (t) =2.5 log (1 



2.5 log 



/ R B (A)$'(A;t) XdX 



Z A 



JR a (X)®' (X(l + z);t)XdX_ 

where A represents the observed filter and B represents 
the rest-frame filter to which we are transforming, Ra 
and Rb are the corresponding observed and rest-frame 
filter response curves, and Za and Z B are the photomet- 
ric zero-points. The K-corrections are therefore applied 
in the following sense: itib = ttla — Kab- The spec- 
tral energy distribution (SED) $' (A;t) at epoch t is ob- 
tained by " color-matching" th e corresponding SED tem- 
plate from IHsiao et alJ (|2007| ). Initially, K-corre ctions 
are computed using unmodified IHsiao et al.l ([2007) SED 
templates and applied to the photometry. Light-curve 
templates are then fit to the N filters for which there 
is optical and NIR photometry yielding N — 1 colors as 
a function of epoch. We then construct a smooth func- 
tion S (A) which, when multiplied by the template SEDs, 
yield synthetic colors equal to the observed colors; i.e., 
we model the template rather than the (more noisy) ob- 
served colors. Finally, these improved SEDs are used 
in Equation [3] to compute the final K-corrections. This 
color-matching simultaneously accounts for the intrinsic 
color variations from SN la to SN la as well as reddening 
c orrections du e to dus t. 

IHsiao et al.1 (|2007l ) demonstrate that while K- 
corrections are determined mainly by broad-band colors, 
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accounting for differences in spectral features is also nec- 
essary. These corrections are based on a larger library 
of spectra, with a greater number of epochs and wave- 
length coverage than previously available. They include 
67 spectra from the CSP (34 of which cover the I-band) , 
obtained at the DuPont telescope, and made available 
to Hsiao ct al. for this purpose. Overall, this new li- 
brary contains many more spectra with red wavelength 
coverage, and the telluric features at A ~ 6880A, 7200A, 
8200A, and 9400A have been identified and removed. In 
the /-band, the Hsiao et al. template includes approxi- 
mately 250 spectra. Determining accurate K-corrections 
in the z-band requires careful attention to the broad Ca 
II triplet P-Cygni absorption feature (8498, 8542 and 
8662 A). Fortunately, for the purposes of the CSP, it is 
only well after maximum light has occurred that the Ca II 
feature changes dramatically ([Hsiao et al.ll2007l see Fig- 
ure 5). Currently, in the spectral region of the i-band, 
125 spectra have been used for the K-correction template. 
We are continuing to acquire additional spectroscopy of 
our low-rcdshift candidates, specifically to improve the 
K-corrcctions function of both epoch and decline 
rate of the SNe la. 

We show in Figure [5] typical cross - band Kjy and Kjj 
corrections based on the lHsiao et al.l (|2007l ) library, plot- 
ted as a function of time since maximum light at B, for 
rcdshifts 0.2, 0.3, 0.45, and 0.6. For comparison we also 
show the K-corrections from lNugent et all (|2002f ). Given 
the currently larger number of available /-band spectra, 
and correction for the presence of telluric features, we 
have adopted the newer K-corrections from Hsaio et al. 
for the purposes of this study. The greatest differences 
occur at early times, 10 days before peak magnitude. 
The corrections are greatest at late times in the SN la 
evolution, as well as very early times for z=0.6. For our 
observations around peak magnitude, the K-corrections 
at these redshifts range from —0.5 to —0.3 magnitudes 
in Y-band and from —1.1 to —0.9 magnitudes in J-band. 
Table[3]coiitains our currently adopted values for the K;y 
and Ki,/ -corrections. In order to estimate the uncertain- 
ties in the K-corrections, we apply the same procedure 
to the SNe la whose spectra were used to generate the 
template SEDs themselves. For those epochs with spec- 
tra, we can compute the K-corrections using the color- 
matched SED and the observed library spectrum, and 
compare the results. For each redshift, we compute syn- 
thetic photometry and K-corrections based on the library 
spectrum and appropriately redshiftcd filter functions, 
simulating what would be observed if this SN la were 
at redshift z. The SED template is then color-matched 
to these synthetic observed colors and the K-correction 
computed again. This is repeated for each library spec- 
trum and the rms difference between the library and 
template K-corrections is computed. Since ESSENCE, 
SNLS, and SDSS-II have three different filter sets, we 
have separately computed the uncertainties in each case. 

In Figure [6] we plot the statistical dispersion in the K- 
corrections as a function of redshift for the filters used 
in ESSENCE, SDSS-II and SNLS. The dispersion is sen- 
sitive to the filter set used, which determines the accu- 
racy with which we can color-correct the template SED. 
As expected, the dispersions are lowest at the redshifts 
where the observed and rest-frame filters overlap and in- 



crease at higher and lower redshifts. The largest uncer- 
tainties arise when there are no filters to anchor the blue 
side of the SED at low redshifts or, alternatively, the red 
side at high redshifts. For comparison, we plot the dis- 
persions with (solid) and without (dashed) including the 
NIR photometry; in most cases, the curves overlap and 
are indistinguishable. In general the statistical uncer- 
tainties in the K-corrections range from ±0.005 to ±0.05 
mag. The largest uncertainty is for dl49 at a redshift 
of 0.34, where the transformation to the Z3-band reaches 
an uncertainty of ±0.1. Similarly, for four SDSS objects 
at z<0.25 (SN 3331, 5549, 7243, and 7512), the uncer- 
tainties also reach ±0.1 mag. For the SNLS objects, the 
statistical uncertainties in the K-corrections are generally 
less than ±0.04 mag. The errors for individual epochs for 
a given object are correlated, and a combined, weighted 
uncertainty is applied to the distance modulus (or peak 
magnitude). These uncertainties are included in the to- 
tal statistical uncertainties computed for each individual 
SN la. We carry these uncertainties in the analysis of 
cosmological parameters described throughout this pa- 
per. Fortunately, progress is continuing to be made in 
improving K-corrections for SNe la. The spectral library 
will continue to grow in the next few years with addi- 
tional spectra from numerous groups, including the low- 
redshift CSP. 

5.2. Extinction Corrections 

Correcting accurately for extinction remains a chal- 
lenge for SN la cosmology. The issue is complicated by 
many factors: 

1) there are at least four separate potential sources of 
dust: a) foreground (Milky Way Galaxy) dust, b) dust 
within the host galaxy, c) dust associated with the cir- 
cumstellar material of the SN progenitor and d) dust in 
the intcrgalactic medium. 

2) there is no a priori knowledge of the dust properties 
in the latter three of these environments and the red- 
dening laws could, in principle, be different in all three 
environments, perhaps being a function of metallicity, at- 
mospheric environment of the SN la, or even evolving as 
a function of time. 

3) even if the extinction law(s) is known perfectly, cor- 
recting for reddening requires knowledge of the intrinsic 
colors (and their dispersion) to derive a reddening curve. 
The observed color of a SN la is determined by the in- 
trinsic SED of the SN la, reddening due to dust, and 
the expansion of the Universe (the K-correction). Hence, 
there is an inherent circularity in the problem, and gen- 
erally some assumptions arc made (e.g., that the intrinsic 
colors of SNe la are known, and that they are known as 
a function of redshift and environment, and/or that the 
reddening law is universal). At present, for SNe la, it is 
still not possible to distinguish unambiguously between 
a different reddening law and differing intrinsic colors. 

The reddening law can be characterized by a ratio 
of total-to-selective absorption, R\, which generally in- 
creases toward shorter wavelengths: R\ = A\/E(B — V), 
where A\ is the total absorption at each wavelength, A. 
On average, the ratio of total-to-selective absorption, R\ 
decreases from 4.9 at U, to 4.1 at B, 3.1 at V and 1. 7 
for the /-band in the Galaxy (e.g.. ICardelli et aLlll989h . 
For a reddening of, say, E(B —V) = 0.02 mag, the cor- 
rections to the rest- frame £/-band magnitude would be 
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~0.10 mag and ~0.03 mag at /. Thus, longcr-wavclcngth 
observations offer a significant advantage in minimizing 
systematic effects due to reddening. 

fn general, it has been concluded that the dust prop- 
erties in other host galaxies appear to be similar to those 
in the Milky Way (see, for example 
Phillips et all [1999 : iKnofj et all 12003 



Riess et al 



Riess et al.l 



1996 



2004 



Jha et al.ll2007t ). To date, no empirical evidence for grey 



dust has be en found fe.g.. lKnop et al.ll2003l : [Riess et al.1 
12004 I2007t ): i.e., larger dust grains with wavelength- 
neutral effects. However, some studies have indicated 
that the reddening law for SNe la is consistent with 
Rv ~ 2.5, lower than the Galactic reddening law (e.g., 
IWang et al.l 1200(1 ). Correcting for SN la reddening has 
generally made use of one or more of the following ap- 
proaches: 1) using a sample of SNe la where the red- 
dening is expected to be negligible in order to define a 
zero-exti nction fiducial sample (e.g., SNe la in e lliptical 
galaxies : lHamuv et al.lll996ctlPhillips et alJl999l ) and/or 
2) making use of a discovery by iLira I ( 19951 ) that after 
~30 days, the ( B — V) colors of SNe la show a very small 
dispersion (e.g. , IPhillips et all 119991 : iPrieto et all 120061 : 
iJha et all l2007f ) , or 3) not correcting directly for red- 
dening, but solving for a general (B — V) color term that 
treats the differences in SN la i ntrinsic color s and redden- 
ing a s indistinguishable (e.g.. iTripp 1119981 lAstier et al.l 
[200l . 

We have dealt with the extinction using two different 
methods. First, we use a "reddening-free" magnitude, as 
described below. For comparison, we also have so l ved fo r 
the reddening explicitly following IPhillips et alj (jl999t ). 
in which an intrinsic color is assumed for a given value 
of Amis and the observed color (after K-corrections and 
Milky- Way reddening corrections are applied) yield the 
reddening due to the host galaxy. We differ from Phillips 
et. al in that we allow for negative reddenings in the 
models. While a negative reddening is not physical, mea- 
suring a negative value for the reddening parameter is 
certainly statistically possible given the uncertainties in 
the measured photometry, and our uncertainty in the in- 
trinsic dispersion of SN la colors. For the whole sample 
of SNe la, the results are statistically unbiased. We list 
the reddenings, corrected for Galactic foreground red- 
dening, for our individual SNe la in Table [5J The mean 
reddening for the low redshift sample is < E(B — V) > 
= 0.06 with a standard deviation of 0.09, consistent with 
that for the high redshift sample, with < E(B — V) > = 
0.05 and a standard deviation of 0.10. 

We proceed to compute a reddenin g-free magnitude, w 
(Ma dori[T982T; iFreedman et al.ll200l[) . defined he re as: 



J BV 



i-R BV (B-V)=i -R BV (B-V) (3) 



where the subscript, 0, refers to intrinsic (unreddened) 
magnitudes and magnitudes without subscripts are ob- 
served magnitudes. The reddening coefficient R BC is 
defined as 



11 BC 



A 



Ra 



E(B -C) R B - Rc 



where R BV 



1.9 and R BI = 0.8 for the case of a stan- 



dard reddening law (with Rv = 3.1), while R BV = 1.1 
and R BI = 0.5 for R v = 2.0 (lCardelli et al.l lT989). The 



advantage of reddening-frec magnitudes is that no knowl- 
edge of either the intrinsic colors of SNe la, nor a sample 
of unreddened SNe la is needed: w l BV is defined such 
that the observed and intrinsic combinations of these 
magnitudes and colors are numerically equivalent. This 
method is in wide use for Cepheid variables. However, 
(just as for other methods), the same color coefficient 
is applied to both the nearby and distant samples. It 
should be n oted that thi s is equ ivalent to using the color 
term, /?, of lAstier et al.l (|2006f ) if the fiducial color of a 
SN la is 0. Indeed, it would simply be a reddening cor- 
rection if all SNe la had zero colors. An advantage of 
this method is in the case where the reddening is solely 
due to dust extinction when this approach will correct 
for it without the need to isolate an unreddened sample. 
Furthermore, the (3 coefficient used by other authors will 
be the reddening coefficient in such a case. 

Given the definition of w BV , and because we are adopt- 
ing a single reddening coefficient, it makes no difference 
to the final results whether we deal with reddening-free 
magnitudes or reddening-corrected magnitudes; the re- 
sults are mathematically equivalent. For both methods, 
we used the high- and low-rcdshift data to determine the 
best value of Ry by minimizing the scatter in the Hub- 
ble diagram, while simultaneously solving for the best-fit 
cosmology. 

5.3. Other Systematic Uncertainties: Evolution, 
Metallicity and Weak Lensing 

In addition to the uncertainties in K-corrections and 
extinction discussed above, there are other potential 
uncertainties on the luminosities of SNe la (e.g., evo- 
lution, metallicity, and weak lensing). For complete- 
ness we briefly summarize the observational situation 
with respect to these effects. It is observed that spi- 
ral galaxies h ost slower decliners, and hence, more l umi- 
nous SNe la (|Hamuv et al.l Il996d Hress et al.lll998[ ). In 
general b righter SNe la occu r in bluer, lower-luminosity 
galaxies (jHamuv et al.|[2000l ). In addition, the scatter 
in the Hubble diagram is observed to be a function of 
morp hol ogical type of t h e hos t galaxy (jSullivan et al.1 
I2003f) . iGallagher et~aTl (|2008f ) have noted a correla- 
tion for nearby E/S0 galaxies such that SNe la in older 
galaxies are fainter than those in galaxies with younger 
global ages, as estimated from stellar population mod- 
els. Further, they find that residuals in the Hubble dia- 
gram correlate with the host-galaxy metallicity Con- 
straints on differing individual SN la properties are 
now coming from detailed comparisons of SN la spec- 
tra dHook et al.ll2005t Igalland et afll2006t iBkmdin et all 
120061 : iRiess et al.l 120071 ). To date, these studies have 
revealed no evidence for significant evolution or metal- 
licity differences. However, there appear to be system- 
atic differences at shorter wavelengths, particularly in 
the re stframe ultraviolet ([Ellis et alJ 120081 : iFolev et all 
|2008f ). Although empirically such effects appear to be 
small relative to the cosmological effect being measured 
(e.g.. lHamuv et aT1ll995l : IRiess et al.1 12004 iKnop et all 
120031 : lAstier et al.1 120061 ) , understanding at what level 
these factors affect the observed properties of SNe la is 
critical to SNe la cosmology. We will return to these 
questions as more data become available and we com- 
plete the CSP. Finally, we note that the net magnifica- 
tion due to weak lensing is not predicted to be significant 
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at redshifts z<0.7 (|Holzlll998D , the redshift interval of the 
CSP sample. 

6. LIGHT CURVES AND DECLINE RATES 

In Figures [7| to [TUl we present optical and Y J-band 
light curves for the 35 SNe la with host galaxy subtrac- 
tions. The Y" J-band data and uncertainties are given in 
Table [3] We also show the optical light-curve data from 
the SNLS, ESSENCE, and SDSS-II surveys. The solid 
lines are template light curves, generated as described 
below. The quality of the light curves in general is quite 
good, and there is an excellent correspondence between 
the optical and the NIR data. For the CSP photometry, 
the average Y-band uncertainties are ±0.03 mag, ±0.06 
mag (±0.10 mag at J), rising to ±0.08 mag (±0.19 mag 
at J) for the redshift intervals 0.1 < z < 0.3; 0.3 < z < 
0.5; and 0.5 < z < 0.7, respectively. 

We make use of our own CSP low-redshift , optical data 
(|Contreras et al.ll2009t iFolatelli et all 120091) for compar- 
ison with our NIR sample at higher redshifts. The 
low-redshift sample is currently comprised of 21 well- 
obs erved SNe la (those l abeled Best Observed in Table 
1 of IFolatelli et all (|2009l )). with redshifts z > 0.01, and 
with E(B — V) < 0. 5 mag. As described in more detail in 
iBurns et al.l (|2009f ), we construct a set of BVgri light- 
curve templates based on this sampl e of 21 CSP SNe la . 
We use a technique similar to that of lPrieto et alJ (|2006l ) 
to generate the light-curve templates shown in Figures [7] 
to 1 101 and determine the Amis decline-rate values and 
time of B maximum using \ 2 minimization. The fits are 
done in flux space simultaneously for the optical and NIR 
data. This technique will be refined as additional data 
are obtained during the low-redshift part of the CSP. As 
shown in Figurc[TJ the j-band light curves usually exhibit 
a second maximum, which can vary in strength from 
event to event. The variations of the i-band template 
are included in the error budget when fitting the z-band 
templates; these become a statistical error when the i- 
band distance moduli are plotted in the Hubble diagram. 
As described in Q5.1[ K-correctio ns were computed by 
first color-matching the SED from IHsiao et alJ (|2007l ) to 
the light curves at each epoch. Independent co-authors 
(CB, GF, and MP) performed a double-blind check of the 
light-curve parameters (Amis), the K-corrections and 
the reddenings. Our derived light-curve parameters are 
given in Tablc[S] As discussed earlier, some of the derived 
E(B — V) values have negative values. Given the mea- 
surement uncertainties and color fluctuations of ±0.06 
mag, however, these negative values are not significant. 

7. A COMPARISON OF PROPERTIES OF THE CSP LOW- 
AND HIGH-REDSHIFT SNE IA 

We compare here the restframc (B — V) and (V — i) 
colors for the CSP nearby and high-redshift samples ana- 
lyzed in this paper. These colors are computed based on 
the peak magnitudes in each band. In Figures [Til and [T2l 
we show histograms for the 21 nearby SNe la (z < 0.1) 
and 35 more distant CSP SNe la with (0.12 < z < 0.70). 
The mean colors and l-a dispersions for our sample are 
(V - i) = -0.63 ± 0.12 mag and (V — i) = -0.67 ± 0.16 
mag at low and high redshifts, respectively. For (B — V). 
the mean colors are (B — V) = 0.04 ± 0.09 mag and 
(B — V) = 0.02 ± 0.10 mag, respectively. The colors are 
consistent to within the uncertainties. However, the very 



red CSP objects SN 2005A and SN 2006X have not been 
included in this analysis since they do not fit our red- 
shift or reddening criteria, as defined above. In Figure 
1131 we show a comparison of the distribution of values 
for Amis for the low- and high-redshift samples. The 
mean values for Amis are 1.11 ± 0.28 and 1.18 ± 0.33 
for the nearby and more distant samples, respectively. 
There is good overlap in the two distributions, although 
the distant sample is more peaked. Given the different 
selection effects for the different samples and the possibil- 
ity of SN la evolution, we do not necessarily expect the 
low- and high-redshift distributions to agree. For this 
small sample, quantitative comparison is limited; how- 
ever, we conclude that the low-redshift sample spans the 
parameter space of high-redshift color and Amis ■ 

8. DISTANCE MODULI, REDDENINGS AND ERRORS 

We now turn to the determination of the distance mod- 
uli, reddenings and the errors associated with these quan- 
tities. We undertake a simultaneous fit for the distance 
modulus as well as three parameters used to define the 
absolute magnitudes. We include both the low- and high- 
z SNe la in this analysis. For each independent choice of 
three filters used to define the reddening-frec magnitude 
(generally BVi), we use the following three parameters to 
define the absolute luminosity of a SN la. We treat these 
three quantities as nuisance parameters in this analysis. 

1) Wq, the absolute reddening-free magnitude of a SN la 
with Am 15 = 1.1 

2) b, the slope of the W-Amisrelation and 

3) Rv, the ratio of total-to-selective absorption. 

The absolute (rcddcning-frcc) magnitude for a SN is 
then Wg V = W + b (Amis - 1-1). To compute the 
distance moduli for the entire sample, we then mea- 
sure an apparent reddening-frec magnitude at maximum 
light, w^y (t max ), for the restframe BVi photometric 
bands. The distance modulus is then, by definition, 
/io = w l BV — W]g V . Distance moduli are listed in Tabled 
(Here, a value of Hq = 72 km/sec/Mpc is adopted. How- 
ever, the value of Hq cancels out for the determination 

of other cosmological parameters in ij9]) 

Usi ng a reduced \ 2 = 1 approach (|Tremaine et alJ 
l2002f ). we then compute x 2 by comparing the measured 
distance modulus, /iq, with a theoretical distance mod- 
ulus ht (C) where C represents the set of cosmological 
parameters (e.g., fl m , &de, w) we consider: 

2 _ V- \H ~ (C, Zj)f 

X ^ a 2 + a 2 K ' 

. cjj -t- a SN 

where &sn is the intrinsic dispersion of SNe la, chosen 
such that reduced-x 2 = 1 (for which we find a value of 
dgjv = 0.09), and a 2 is the total variance for observation 
j and is given by 

a 2 = a 2 (i max ) + R l B V cr 2 {B max - V max ) + b 2 a 2 (Ami 5 ) 
2i?^jy(r(i mQ2 ;, B max V max ) — 1bo(i max , Amis) 

in ax 

-V max ,Am 15 ) (5) 

where cr 2 (x) is the variance in parameter x and a(x,y) 
is the covariance between parameters x and y. A pecu- 
liar velocity term of ±300 km/sec is also included. By 



10 



Frccdman ct al. 



minimizing this \ 2 (i.e., the residuals in the Hubble dia- 
gram), we simultaneously determine the cosmology and 
the three nuisance parameters Wq, b, and Ry. In this 
sense, our method is simil a r con c eptually to that de - 
scribed by [Tripp fc Bra nd] (|1999ft : lAstier et all (|2006ft : 
IConlev et alJ (|2008f ) . We defer a discussion of the re- 
sults for cosmology until SjHl and discuss first the results 
for the decline-rate and reddening parameters from this 
method. 

Using the i-band data at maximum, corrected for red- 
dening based on the BV photometry, the results for W B y 
are W = -18.45 ±0.05 (statistical) ±0.01 (systematic), 
6 = 0.38 ± 0.08 (statistical) ±0.01 (systematic), and 
R v = 1.74±0.27 (statistical) ±0.1 (systematic). 26 Since 
we have data for three filters, we can independently cor- 
rect for reddening using different combinations of band- 
passes (e.g., B — V, V — i, or B — i). However, we prefer 
to solve for the reddening using simultaneous data from 
a single photometric/telescope system (B — V), which 
is available for all of the SNe la, rather than a hybrid 
optical-NIR combination. If instead we use the -B-band 
data and correct for reddening based on the BV photom- 
etry W§ v , we find Wo = -19.11 ±0.07, b = 0.65 ±0.13, 
and Ry = 1.66 ± 0.27. Encouragingly, Ry is consis- 
tent for the two filter combinations. The final results for 
the cosmology are consistent to within the uncertainties, 
with all combinations of bandpasses. 

A value of Ry = 1.74 is significantly lower than 
a Milky- Way value of Ry = 3.1. A lower value of 
Rv is also found by iTripp I (fl998h. ITripp fe Branch! 
(11999ft. andlWang et all (120061). IConlev et all (|2008f ) and 



lAstier et al j (|2006ft would find an even lower value than 
ours, consistent with Ry (= Rb — 1) ~ 1 (their j3 = 2 
would correspond to Rb), if their color term was associ- 
ated with the extinction law. 

As a check on our reddening corrections, two indepen- 
dent analysis methods were used by two of this paper's 
co-authors (CB and GF) to determine the nuisance pa- 
rameters: the reddening-free meth od outlined above, an d 
the extinction method outlined in lPhillips et alJ (|1999ft . 
We initially discovered a large discrepancy between the 
two derived values of Ry. The reason for this discrep- 
ancy can be understood from equations 2] and [5j The 
denominator of \ 2 includes terms that depend on both 
Ry and b. Depending on the magnitude of the variances 
and covariances, these terms have leverage on the final 
solution. We have found that a significant issue affect- 
ing the derived reddening law is what is assumed for the 
variance in SNe la colors, a 2 (B — V), with the value of 
Ry increasing with increasing a 1 (B — V) (sec Appendix 
|A")) . The two analysis methods agree well if consistent 
errors are adopted. It should therefore be stressed that 
accurate estimates of the variances in the data are crucial 
for the determinatio n of Ry. As a furth er test of these 
different approaches. iBurns et al. I (|2009ft have begun to 
explore this issue with an unbiased estimator using a 
Monte Carlo Markov Chain (MCMC). This work goes 
beyond the scope of the present paper, but we note that 
the preliminary results agree well with those presented 
he re, with Ry = 1.8. Our current analysis and that 
of iFolatelli etaLl (120091) are consistent with a picture in 

26 Unless otherwise noted, all reported uncertainties are statis- 
tical and correspond to 1-cr errors (68% confidence). 



which, in addition to corrections for Milky- Way- like dust, 
there is an intrinsic dispersion in the colors of SNe la that 
is correlated with luminosity, but not the decline rate. 

Before concluding this section, it is worth emphasizing 
that disentangling a possible additional intrinsic color 
term from extinction leaves open a challenge to SN la 
cosmology. Potential differences in cither the intrinsic 
or observed sample populations of SNe la with redshift 
and/or environment, or differences in the reddening law 
may limit the ultimate accuracy for cosmology, and for 
reaching systematic uncertainties significantly below the 
~ ±10% level. Having highlighted the remaining chal- 
lenges, we note, however, that the level of the effect we 
are discussing amounts to less than ±10%, and reach- 
ing a systematic uncertainty at this level for SNe la is 
already excellent progress. Few other methods in use 
for dark energy measurements currently offer this kind 
of accuracy. We now turn to a discussion of the Hubble 
diagram, and then to implications of different values of 
Ry on the cosmology. 

9. HUBBLE DIAGRAM AND COSMOLOGICAL RESULTS 

9.1. The First I-band Hubble Diagram for z > 0.1 

In Figure [2J we show the i and i?-band Hubble dia- 
grams, uncorrected for reddening. The residuals about 
the best-fit to these data are shown in the panels be- 
low. We can immediately see the advantage of observ- 
ing SNe la at red wavelengths relative to the optical. 
The rms dispersions in these plots amount to ±0.17 and 
±0.29 mag, respectively 

In Figure 1151 we present the first i-band Hubble dia- 
gram for the CSP sample of 35 SNe la based on Mag- 
ellan data (solid blue squares) using the reddening-free 
magnitude technique described in ^ adopting a value 
of Ry = 1.74 ± 0.27. The error bars shown in the lower 
panel represent the formal 1-cr uncertainties in the dis- 
tance modulus and include the combined errors in the 
photometry, the dispersion in the light-curve templates, 
the estimated dispersion in K-corrections, the error in 
the foreground Galactic reddening, and the co-variances 
between the light-curve parameters. We also include a 
peculiar- velocity component of ±300 km/sec. The top 
curve corresponds to an Q m = 0.27 , SIde = 0.73 cos- 
mology. For comparison, the small black sq uares indi- 
cate distance moduli from lAstier et a l. (2006) . The cur- 
rent C SP sample of 21 nearby SNe la from lFolatelli et aTl 
(|2009ft (solid red circles) is also shown in Figure [15] The 
low-redshift sample is restricted to nearby SNe la in the 
Hubble flow, having redshifts greater than z > 0.010, 
so that the scatter due to random peculiar velocities is 
minimized. (The results remain consistent, to within the 
quoted uncertainties, if we further restrict the sample to 
z > 0.015 or z > 0.02.) 

The immediate conclusion we can draw from this Hub- 
ble diagram is that the CSP data alone provide indepen- 
dent evidence for a (standard) cosmological model with 
O m = 0.3 and Q,de = 0.7. The CSP constraints can 
be further improved by combining them with other in- 
dependent measure ments, for example, ba ryonic acous- 
tic oscillations fe.g.. lEisenstein et al.ll206"5ft . A weighted 
fit to our 35 data points in combination with baryon 
acoustic oscillations (assuming w = —1) yields the so- 
lution: n m = 0.27 ± 0.02 (statistical), n DE = 0.76 ± 
0.13 (statistical) ± 0.09 (systematic). The systematic 
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uncertainties for SNe la do not have much impact on the 
determination of f2 m since this parameter is determined 
largely by the matter power spectrum. We quote only 
statistical uncertainties for D. m here. The statistical un- 
certainties are determined by marginalizing over all other 
parameters and fitting the ID probability distribution to 
a Gaussian. The systematic uncertainties are treated be- 
low in SgO] 

Based on the fit to the Hubble diagram above, we 
show, in Figure 1161 our error ellipses in the fl m -fluE 
plane. Consistent with previous SN la studies, we find 
that based on the SN la data alone (and the assumption 
that w = —1), a value of £Ide > is required at greater 
than the 99% confidence level. 

As an alternative to making assumptions about if, 
we can use the CSP data to calculate w under the as- 
sumption of flatness (flk = 0). Here we again com- 
bine the CSP results with independe nt measurements 
of ba ryonic acoustic oscillations (e.g., lEisenstein et alJ 
120051 ) , as shown in Figure [TTl Assuming a flat cosmol- 
ogy, these joint constraints yield a value of f2 m = 0.27 ± 
0.03 (statistical) and w = —1.05 ± 0.13 (statistical) ± 
0.09 (systematic). For the purposes of this calcula- 
tion, we are assuming that w is a constant (i.e., w a = 
0). These results are in excellent agreement with other 
joint constr aints from SN la studies and baryon acoustic 
oscillations (lAstier et all 120061 : IWood-Vasey et all 120071 
iRiess et all I2007D . which also yield values of w = — 1 
and Q m = 0.3, to within the quoted measurement uncer- 
tainties. We have also combined our CSP results with 
the two-dimensional pr obability contours fro m the 2dF 
galaxy redshift survey (|Hawkins et all l2003f ) . We find 
tt m = 0.27 ± 0.09 (statistical) and w = -1.03 ± 0.23 
(statistical), in good agreement with the SDSS data, but 
with larger uncertainties on the value of wq. 

In Figure [TBI we show a comparison of the residuals in 
the Hubble diagram for the high-z data (z > 0.10) for 
the i- and i?-bands relative to the best-fit flat, constant 
w cosmological model, again for a value of Ry = 1.74. 
The rms scatter in the i-band Hubble diagram amounts 
to ±0.13 mag about the fit ; the B-band scatter amounts 
to ±0.15 mag. 27 The resulting cosmological parame- 
ters are in good agreement, with differences in w and 
O m amounting to only 0.03 and 0.01, respectively, well 
within the measurement uncertainties. The z-band mea- 
surements have smaller systematic uncertainties due to 
smaller reddening corrections, and also have somewhat 
smaller observed scatter. 

As discussed in Sj6l we find that the smallest scatter in 
the Hubble diagram at low (as well as high) redshift is 
found for a value of the (color coefficient or) ratio of total- 
to-selective absorption Ry = 1.74 ± 0.27. The standard 
Galactic reddening law is Ry =3.1. In Table [3 we 
show the fits to Vt m and £Ide and the rms scatter about 
the cosmological fit for two values of Ry, 1.74 and 3.1, 
for both the i and B-band data. For an adopted value 
of Ry = 1.74, the resulting values of fi m and w are 
in excellent agreement between both bands. However, 
adopting instead a value Ry = 3.1 results in a value of 
w = -1.20 ± 0.13 for the i data and w = -1.24 ± 0.16 

27 The rms values excludes SN 04D2an at a redshift of z = 0.62, 
which has no reliable rest-frame V-band observation and therefore 
no color correction is possible using only the optical data. 



for the B data. The solutions for these different adopted 
values of the reddening are still mutually consistent, but 
the scatter in the Hubble diagram increases by about 
50% for Ry = 3.1. Given the increase in the Hubble 
diagram scatter at larger redshifts for values of Ry = 3.1, 
and the excellent fit to the low redshift CSP sample for 
Ry = 1.74, we conclude that currently the CSP data are 
consistent with a lower value of a color coefficient and/or 
reddening law. 

It is interesting to ask what happens if reddening ef- 
fects are ignored altogether. In this case we find (based 
on the data shown in Figure IT4")) . that for i, fl m = 0.29 
± 0.03 (statistical) and w = -0.90 ± 0.14 (statistical), 
consistent to within the systematic uncertainties with the 
results from the reddening-frec and reddening-corrected 
data discussed above. For the i?-band data, however, f2 m 
= 0.31 ± 0.03 (statistical), w = -0.70 ± 0.21 (statisti- 
cal), and the results for the equation of state are much 
more sensitive to the treatment of reddening. These re- 
sults illustrate quantitatively the advantage of the i band 
in minimizing systematic effects for SNe la cosmology. 
We have also shown that the i-band can be calibrated as 
well as the B-band; in fact, there is a further advantage 
to the i-band, where the luminosity-Amis relation has a 
smaller slope. 

It is beyond the scope of this paper to ascertain unam- 
biguously the reddening properties for SNe la; however, 
we make the following observations. The simplest as- 
sumption, the application of a standard reddening law 
coefficient, Ry = 3.1 alone, docs not provide the best fit 
to either the low redshift or high redshift CSP data as de- 
fined by reduced scatter in the Hubble diagram. As noted 
previously, this in and of itself does not imply a differ- 
ent reddening law within the host galaxies. For example, 
it could be indicating an intrinsic color term or perhaps 
an unusual kind of reddening resulting from dust in the 
vicinity of the SN la. In fact, it would be a somewhat 
unusual circumstance for the dust within the interstellar 
medium of all SN la host galaxies to have different prop- 
erties from the dust in nearby galaxies; e.g., LMC, SMC 
and the Milky Way, which are the only galaxies for which 
the value of Ry can be measured to high accur acy and 
wher e on average it is found to b e equal to 3.1 (|Draind 
I2003D . iNobili and Goobar I (|2008l) have determined the 
average reddening l aw for a sampl e of 80 SNe la and find 
that it agrees with ICardelli et all (|1989l ) for R v 1.75 ± 
0.27. It has also recently been shown that normal dust 
which is distributed locally around th e SN can produ ce a 
reddening law with lower value of Ry (|Goobar 1120081 ). In 
other words, this abnormally low value of Ry may sim- 
ply be due to the geometry of local dust around the SN, 
which dominates any extinction produced by the host 
galaxy. 

We plan a detailed investigation of the reddening law 
for SNe la based on our low redshift CSP sample, which 
was designed to provide independent determinations of 
reddening at multiple wavelengths. For the present dis- 
cussion, we adopt a general color coefficient term con- 
sistent with a reddening law (or color term), Ry = 
1.74 ±0.27. 

9.2. Dark Energy and Other Parameterizations of the 
Expansion 
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How do we best compare the observational SN la data 
with cosmological models? While a simple question, this 
issue is non-trivial. As discussed in iJTJ the expansion of 
the Universe can be parameterized in terms of the equa- 
tion of state parameter, w. Lacking a physical explana- 
tion for the dark energy, this is a helpful parameteriza- 
tion, but it is still just that - a parameterization. It is 
important to keep in mind what the actual observables 
are: observations yield magnitudes and redshifts - not 
accelerations, equations of state or second derivatives of 
scale factors (the deceleration parameter) . In the case of 
a flat universe, the luminosity distance is given by: 

d L (z) = r(z)(l + z) = £ j^—, (6) 

where r(z) is the comoving distance. Because the lumi- 
nosity distance relates (inversely) to the integral of H(z), 
it does not provide information on how the expansion 
rate may have changed from higher redshifts to today. A 
model is required to infer this evolution. 

We now turn to a discussion of the equation-of-state 
parameter, w, which is commonly used in current cos- 
mological models. A number of ways to pa rameterize w 
have been suggested in the literature fe.g.. lLinderll2006l ; 
lAlbrecht et al.ll2006f ). The functional forms for the evo- 
lution in w have included simple terms that are either 
linear functions of rcdshift, z (or scale factor, a), but 
range also to more exotic functional forms. It must be 
emphasized that while SN la data alone provide a very 
strong case for a non-zero value for Qde (i-c, an addi- 
tional component to matter and radiation), SN la data 
alone do not provide a sensitive means of constraining 
either the value of w at the current epoch, or its time 
evolution. 

As discussed in Q9.1\ under the assumption of flatness 
(Qk = 0), SNe la can provide a strong constraint on 
w when combined with another experiment like baryon 
acoustic oscillations. Given current evidence, flatness is 
not an unreasonable assumption; for example, we note 
that either recent Hq or SNe la results, in combination 
with the WMAP ob servations, yields fi fc = -0.01±0.01 
(|Spergel et al.ll200"7t ). However, additional methods (e.g., 
baryon acoustic oscillations, cluster growth, weak lens- 
ing) must be combined with the SN la data to provide 
meaningful constraints on w. At the current time, the 
accuracy required is not yet sufficient for any of these 
methods alone, and is the goal of future studies. 

The cosmological parameterization (f2 m , Qde, w) rests 
on assumptions about the matter/energy content of the 
Universe and the functional form for the dark energy 
equation of state. Here we also investigate another pa- 
rameterization that is independent of this theoretical 
framework and involves a purely kinematic model de- 
scribed by the 3 parameters: (go,jo,^fc)- Here go = 
— aa~ 2 a is the cosmic deceleration; jo =a a~ 3 a 2 is the 
third derivative of the scale factor, the so-called cosmic 
jerk; and flk is the curvature parameter. The only as- 
sumption that enters into this parameterization is that 
the Robertson- Walker metric accurately describes the ge- 
ometry of the Universe. It is therefore free of assump- 
tions about the energy content of the Universe or even 
that the Einstein field equations are the correct descrip- 
tion of gravity; it is a purely kinematic model. Although 



currently somewhat out of fashion, it still remains the 
closest parameterization of the data, since it is acceler- 
ation that is actually being measured. In the standard 
cosmological model (Cl m , Ude, w) = (0.3, . 7, — 1) , j = 1, 
and q = -0.67. As noted bv iBlandforl (|2004f ). lack- 
ing an understanding of the dynamics of the Universe, 
a purely kinematic description remains a well-motivated 
family of models to explore. 
The scale factor is expanded as a Taylor series: 

a{ tt\ o) ~ l+H (t-t )-^q H 2 (t-t ) 2 + y o H^t-t f+. . . 

(7) 

from which one can then d erive a luminosity distance. 
ICaldwell fc KamionkowsO (|2004h provide a convenient 
expansion of the luminosity distance to third order in z 
and show that to within the precision of this truncation 
of the Taylor series, the cosmic jerk and curvature can 
be combined into one parameter (jk = j + Qk)' 

, , * cz ( 1 . 

L (z) "i5o \ + 2 ( 90)Z_ 

i (l-g - Mo+Jo + ilk) * 2 + C* 3 )} • 

Using the definitions of these parameters and the Fried- 
mann equation, one can derive the transformation equa- 
tions between the two parameterizations: 

g = i (ft m + n DE (1 + 3io)) 

jo = n rn + ^(2 + 9w(l + w)) (8) 

Ofc = Q m + Q.DE — 1 

In Figure I19[ we show the sum of the jerk and cur- 
vature (j + Qk) parameters as a function of go- A con- 
stant value of the jerk is assumed. The grey shading 
indicates the region where the luminosity distance ex- 
pansion is vali d. 28 The blue contours represent the CSP 
data using the ICaldwell fc Kamio nkowskil (|2004f ) param- 
eterization. The black contours, representing the baryon 
acoustic oscillation data, were generated using the trans- 
formation equations above. The joint constraints from 
the CSP and baryon acoustic oscillation data yield a 
value of jk = 1.18 ±0.44 (statistical) ±0.27 (systematic) 
and q = —0.67 ± 0.13 (statistical) ± 0.09 (systematic) 
at the 95% confidence level. Since jerk and curvature 
are combined, no assumption of flatness is required. We 
note that to within the uncertainties, these parameters 
are consistent with the standard dynamical model with 
w = — 1. With future, larger data sets, the discrimina- 
tion amongst competing models will be sharpened and 
this type of kinematical prescription will offer a valuable 
independent test of the current standard cosmology. 

9.3. Systematic Uncertainties 

To quantify the effects of systematic errors, we model 
the effects of known and potential systematic errors for 
each of the five cosmological parameters discussed in this 
paper (w, Cl m , CIde, qo, an( l Jo)- We present in Table 
[8] the main sources of uncertainty that could contribute 
to a systematic error, including magnitude and color off- 
sets between the low- and high-redshift samples, errors 

28 Specifically, the error incurred by not including the 4th order 
term in the Taylor series expansion results in a distance modulus 
error of ±0.3 mag at a redshift of 0.7. 
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in the color terms used to transform the instrumental 
magnitudes to our natural system, and biases due to the 
method of ^-minimization. We have aimed to give a 
conservative estimate of the uncertainties. The first col- 
umn describes the potential source of error and the sec- 
ond column provides a bound on its magnitude. For each 
of the cosmological parameters, we tabulate the rate of 
change of the parameter with respect to the systematic 
(columns 3, 5, 7, 9, 11, 13, and 15) and the resulting 
error on the parameter (columns 4, 6, 8, 10, 12, 14, and 
16). The main sources of potential systematic errors and 
how we simulate their effects are discussed in more detail 
in Appendix lAl 

The final systematic errors adopted are obtained by 
summing in quadrature the contributions listed in Table 
[5J From this table, we can see that the current systematic 
total uncertainty in our measurement of w is ±9%. A key 
element for future SNe la studies is improving the abso- 
lute calibration of photometric standards (particularly in 
the era of the Joint Dark Energy Mission, JDEM), ensur- 
ing that the calibration minimizes color uncertainties as 
optical and NIR measurements are compared. Decreas- 
ing the uncertainty due to reddening is another critical 
component of minimizing the overall systematic errors. 

Finally, we summarize in Table [9] the values of cos- 
mological parameters from this paper calculated under 
different sets of assumptions, as described in tj9.ll along 
with both their statistical and systematic uncertainties. 
When the analysis of our total sample of 75 high-redshift 
SNe la is complete, our statistical and systematic uncer- 
tainties will be comparable. 

10. SUMMARY AND FUTURE MEASUREMENTS 

We have used ground-based NIR measurements of 
SNe la to yield an independent Hubble diagram based 
on rest-frame z-band data. Reddening effects are a pri- 
ori lower in the red than the ultraviolet-blue-visual, and 
we find that i-band photometry is an effective tool for 
mimimizing systematic effects for SNe la. Our new ob- 
servations of 21 SNe la at 0.01 < z < 0.08 and 35 SNe la 
at 0.12 < z < 0.70, yield the following results: 

1) These first CSP data provide independent evidence for 
an accelerating universe. In the context of a cosmological 
model including a component of dark energy, £Ide > 
at signficantly greater than the 99% confidence level. 

2) Joint SN la plus baryon acoustic oscillation constraints 
yield w = —1.05 ± 0.13 (statistical) ± 0.09 (systematic) 
and O m = 0.27 ± 0.03 (statistical) . When the analysis 
of our total sample of 75 high-redshift SNe la is com- 
plete, our statistical and systematic uncertainties will be 
comparable, and at the ±10% level. 

3) A purely kinematic solution, with no assumptions 
about the matter and energy content of the Uni- 
verse, yields values of jk = 1.18 ± 0.44 (statistical) ± 
0.27 (systematic) and q = — 0.67± 0.13 (statistical) ± 
0.09 (systematic) for the cosmic jerk and the decelera- 
tion parameter, respectively. These results are consistent 
with an acceleration of the expansion of the Universe, 
and with the current standard model of cosmology. 

4) The current sample of SN la photometry is inconsis- 
tent with the application of a standard Milky Way red- 
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dening law alone, suggesting either that intrinsic color ef- 
fects dominate the standard reddening corrections, that 
the SNe la are being reddened by dust with different 
properties than that in the Milky Way, or perhaps there 
is circumstcllar dust about the SN la. Future SNe la 
studies will need to disentangle these effect to decrease 
the systematic errors for SN la cosmology. 

Beyond a redshift of 0.7, the z-band (at restframe 
7600A) is shifted beyond the \.2\ira J-band in the 
NIR, rendering i-band measurements impossible from 
the ground. Even before this limit is reached, observa- 
tions from space are desirable in order to eliminate the 
bright terrestrial sky background. The CSP data illus- 
trate quantitatively the utility of longer-wavelength data 
in minimizing systematic uncertainties for SNe la cosmol- 
ogy, of relevance for future planned space missions such 
as the NASA/DOE Joint Dark Energy Mission (JDEM). 
We suggest that a combination of ground, HST, and fu- 
ture space observations are needed to measure accurately 
the i-band Hubble diagram and constrain the values of 
w and w a to the highest possible accuracy. 

Finally, we note that while SNe la currently provide the 
most compelling evidence for the acceleration of the Uni- 
verse, ultimately a combination of different techniques 
will be required to measure not only w, but also its time 
evolution. Existing degeneracies are such that accurate 
constraints can be obtained only in the combination of 
several techniques, or explicitly making the assumption 
that the Universe is flat (VLtot = !)■ Future studies 
of baryon acoustic oscillations, weak lensing, clusters of 
galaxies, and SNe la will yield further complementary 
and independent estimates of w and w a . SNe la will re- 
main a valuable component of future investigations be- 
cause they naturally cover the redshift range (0.1 to 0.6) 
where dark energy is measured to be dominant, and cos- 
mic variance is not a major issue for them. 
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APPENDIX 
SYSTEMATIC ERROR BUDGET 

In this appendix, we outline our method for estimating the effects of systematic errors in the determination of various 
parameters. We have categorized the various sources of error based on their functional form. We then compute the 
sensitivity of the cosmological parameters to these functions through simulations. Finally, we compute the magnitude 
of each effect by summing the systematics for each category in quadrature. Table [5] summarizes the sensitivity of each 
parameter to each systematic. We also illustrate the effects of these systematic effects on the calculated parameters 
in graphical form. 

Magnitude Offset Between Low- and High-z 

The low- and high-z observations are obtained on different telescopes, and therefore are calibrated independently, 
which can potentially introduce a simple magnitude offset between the zero-points of the low- and high-z data. We 
model this effect by simply adding fixed values, 5 m , to the high-z peak magnitudes and re- compute the best-fit 
cosmology. 

The results of this simulation are shown in Figure [2"01 and show that to a good approximation, the systematic effect 
of a magnitude offset on w is a linear function of 5 m with slope dw/dd m = —2.69 in the i-band and dw/d6 m = —2.66 
in the i?-band. Furthermore, there is no significant effect on the 3 nuisance parameters. 

To estimate the value of S m , we consider the following possible sources of such an offset: 1) error in the Y and J 
zero-points, 2) error in the i-band zero-point, and 3) errors in the extinction coefficients. We describe each in turn. 

The contribution to the systematic uncertainty in S m has a different origin for the Y and J filte rs. As described 
in 14. H we have used Kur ucz models to exten d the energy distributions of the IPersson et al.1 (|1998l) standard s from 
J (1.25/i) to Y (1.035/i) (jHamuv et al.l 12006). That this was viable was confirmed by IContreras et al.l ()2009f ) , who 
showed that the zero-color AO star Feige 16 does in fact have Y — J — —0.009 ± 0.016, confirming this calibration to 
an accuracy of ±0.01 mag. Conservatively, we adopt an error in the zero-point for the V-band of ±0.025. 

In the case of J, the main problem is the changing amount of water vapor above the telescope. An increase manifests 
itself as a systematic narrowing of the passband. This may or may not show up in standard star measurements, 
depending on when during the night they were mea sure d. We adopt an err or in the zero-point for J-band of ±0.015. 

The common ancestor to the IPersson et al.1 ()1998f ) and lSmith et al.l (|2002l ) systems is Vega. We therefore must deter- 
mine the zero-point of the i-band relative to Vega instead of BD + 17°4708 in order to compute consistent zero-points. 
We have i nvestiga t ed th e effe cts of a systematic error in the SED of Vega by comparing synthetic i-band photometry 
using the iBohlinl l)2007f) and lBohlin fc Gillilandl (|2004h SEDs. The difference is 0.01 mag. We further investigated 
uncertainties in our i-band filter transmission function. In constructing all our CSP filters, s hifts in wavelength ne eded 
to be applied in order to have the synthetic color terms match the observed color terms (|Contreras et al.ll2009t ). A 
typical uncertainty in these shifts is approximately lOA, which is equivalent to an error of ±0.001 mag. We therefore 
estimate a zero -point error of ±0.01 for t he i-band. 

According to lHamuv et al.l (|2006f ) and IContreras et alJ (|2009h , the dispersion in the extinction coefficent for the i- 
band is ±0.03. The median airmass of our SNc la observations in the i-band was 1.315 and so we estimate a systematic 
error of ±0.009 mag due to dispersion in the extinction coefficient. Similar errors are expected for the Y- and J-band 
filters, though the median airmass for our high-redshift observations is typically lower (1.24). We therefore estimate a 
total error of ±0.01 mag. 

Given the external check with the standard star Feige 16 (jContreras et al.l I2009T ) . we believe our calibration to 
be robust. Nevertheless, wc allow for a conservative systematic magnitude offset between low- and high-redshift of 
S m = ±0.025 mag. 

Color Offset Between Low- and High-z 

This systematic is analogous to the magnitude offset and can occur for the same reasons; however, it propagates 
through the analysis differently as it is multiplied by the reddening coefficient. We proceed in the same manner as 
before, adding a color offset 5 C to the (B — V) colors at high-z and compute the best-fit cosmology. As in the previous 
case, the effect on the cosmological parameters is linear (see Figure [2Tj) . However, the effect on the B-band is more 
than twice that in the i-band (dw/dS c = 2.53 and dw/dS c = 7.53, respectively). This is simply due to the fact that 
the reddening coefficient is twice as large at B relative to the «-band. An advantage is that the colors we use in our 
analysis are all constructed from the same optical photometric systems and avoid the NIR and i-band. We estimate 
the total possible color offset to be S c = ±0.02. 

Color Gradient 

Color terms are used to transform the magnitudes of the standard stars to the natural CSP system. These color 
terms are empirically determined and therefore suffer from measurement uncertainty. An error in the color terms 
could potentially introduce a color gradient in the data: color errors that correlate with the instrumental color of the 
SNe l a. Figure [23| shows t he se nsitivity functions. As expected, the impact in the B-band is larger than in the i-band. 
From lHamuv et al.1 (|2006l ) and IContreras et all (|2009l ). we estimate 8(CT) = 0.015, which is the typical 1 -a error m 
the color terms. 
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Estimation of Uncertainties 

The me thod of ^-m inimization suffers from a well-known bias when there is measurement error in the independent 
variables (Kelly 2007). The root of the problem lies in the fact that the numerator of \ 2 ( scc cquation[4]) is sensitive to 
the model parameters only, whereas the denominator is sensitive to the nuisance parameters as well as the variances 
(see equation[S|). As such, the best-fit solution is a function not only of the data, but the variances as well, complicating 
the computation of the overall systematic uncertainties. 

We have investigated the effect of adding extra variance to both the colors a\ and Am.15. Figure l2"2l shows the 
strong dependence of Ry on a c . The i-band data are about a factor of two less sensitive to changes in Ry than the 
-B-band data. As a result, the bias due to the denominator of \ 2 is larger in the z-band. If we were to include more 
heavily reddened objects (for instance SN 2005A and SN 2006X), then this effect would disappear entirely. The effects 
of extra variance in Amis are not significant and are therefore not shown, though they are included in Table [8] for 
completeness. 
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Fig. 1. — A representative sample of four i-band lightcurves from the low-z CSP sample chosen to illustrate the behavior of the secondary 
maximum. For clarity, the lightcurves have been normalized to their peak magnitudes and then offset by one magnitude from each another. 
The SNe la have values of Amis ranging from 1.12 to 1.81. The fast-declining SN la with Amis = 1-81 (bottom light curve) shows no 
secondary maximum. 




6000 7000 8000 9000 

Wavelength (Angstroms) 

Fig. 2. — Filter response functions for the Y (blue) and J (red) bands blucshiftcd by 0.35 and 0.63, respectively. Rest-frame i-band is 
slotted in black. T he absorption feature at A ~ 7600A is telluric 02- Crossing the entire figure, the SED of a typical SN la at maximum 
Hsiao ct al. 2007) is also shown in black. 
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Fig. 3. — Redshift distribution of objects observed as part of the CSP in the NIR on Magellan. In general, the higher-redshift objects 
are found by the SNLS, and those at lower redshifts primarily by the SDSS-II. We have observed 14 objects discovered by ESSENCE, 37 
SNLS objects and 24 SDSS objects. 
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Fig. 4. — Sample of SNe la, post-SNe la and difference frames for three SNe la: SN 3241 (discovered by SDSS-I1, top panel), 04Dlrh 
(discovered by SNLS, middle panel), and 05D2bt (discovered by SNLS, bottom panel). SN 3241 (z = 0.25) was observed in the Y-band, 
04Dlrh (z = 0.435) was observed in both Y and J (J is shown), and 05D2bt (z = 0.679) was observed in the J-band. The scale of these 
images is 200 pixels = 25 arc-seconds on a side. 




Fig. 5. — Cross-band K-corrections (Kij and K;y), based on SED templates of lHsiao et al.l {2007]) and Nugent et al. (2002), shown as 
a function of epoch for redshifts z = 0.20, 0.30, 0.45 and 0.60. The solid lines are the Hsiao et al. corrections and the dashed lines are 
Nugent et al. The largest differences occur at about 10 days before maximum. Given the significantly larger template library available to 
Hsiao et al., and the corrections for the telluric features, we have adopted the Hsiao et al. corrections. 
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Fig. 6. — Estimated dispersion in the K-corrections as a function of redshift for different photometric filter sets. Solid lines represent 
the dispersion when full optical and NIR filter sets are used to color-correct the template SED, whereas dashed lines show the dispersions 
when no NIR data are included. The dispersion in the Kftg for ESSENCE increases dramatically at low redshift because Rn m is redder 
than rest-frame B and there are therefore no observations to anchor the SED template on the blue side to which we are transforming. Our 
lowest redshift ESSENCE object in this analysis is at z = 0.34. At the high-redshift end, the uncertainty for Kjy is higher since is 
shifted to the blue side of rest i, and there is no anchor on the red side. Adding the NIR photometry helps in the color-matching and 
decreases the dispersion. There is little effect on the SNLS dispersions because our CSP sample includes only objects with redshifts less 
than 0.8 and the z m band serves as a red anchor over this entire range. 
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Fig. 7. — Optical and YJ-band light curves for SNe la discovered as part of the SNLS, ESSENCE and SDSS-II projects, 
followed up using PANIC on Magellan. Light curves and templates are shown in the observ er (not rest) frame. The template 
BVri light curves (solid lines) are based on the low-redshift CSP data (jFolatelli et al.|[2009l ). and are expanded by (1 + z) and 
K-corrected to fit the observed light curves. The curves are labelled BVri to indicate which restframe template has been used 
to fit the observed data. An offset of 1 mag is applied so that the curves do not intersect. In some cases, both the observed 
Y and J data are de-redshifted to the i band. The maximum-light magnitudes for each filter are presented in Table [S] The 
red and blue stars correspond to the CSP's J- and Y-band data, respectively. The black circles correspond to g m (SNLS) or g s 
(SDSSII); the black squares correspond to r m (SNLS), i s (SDSSII), or R4 m (ESSENCE); the black triangles correspond to i m 
(SNLS), is (SDSSII), or I 4m (ESSENCE); and the black diamonds correspond to z m (SNLS). 
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Fig. 8. — Same as Figure [7] 
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Fig. 9. — Same as Figure [7] 
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Fig. 10. — Same as Figure [7] 
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Fig. 11. — Histog rams of the (V — i) color distributions at low and high redshifts for the SNe la labeled best observed in Table 1 of 
IFolatelli et al. (2009)), with redshifts z > 0.01, and with E(B — V) < 0.5 mag. Note that the (V — i) colors on our natural system are all 
negative, and bluer than (V — I). The object SN 04D2an is not included in this plot as it has no rest- frame V observation. 
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Fig. 14. — Top panel: i and B-band Hubble diagrams for 21 low-redshift and 35 high-redshift SNc la from the CSP, uncorrected for 
reddening. Bottom panel: The residuals about the best-fit to these data. The values for rms scatter about the best fit to these data are 
labelled. The rms value in brackets excludes the most discrepant (highly reddened) SNLS 05Dlhn. 
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Fig. 15. — j-band rest-frame Hubble diagram for a total of 5 6 CSP SNe la, 35 SNe la from the Magellan CSP sample (blu e squa res), 
and 21 low-redshift data (red solid circles) from Folatclli ct al. (2009). For comparison, distances determined by [Asticr ct al. (200j]|) are 
shown as black squares, but are not included in the fits. Error bars shown are 1-ct. A value of Ho = 72 km s — 1 Mpc" 1 has been adopted 
for the plot. The solid (red), dot-dashed (black), and dashed (blue) lines represent Q m = 0.3, Orj_g = 0.7 ; Q m = 0.3, £Ide = I an d 
= 1 models, respectively. The data are consistent with the standard (accelerating) cosmological model. To minimize the effects of peculiar 
velocities, the fit to the low-redshift sample is restricted to z>0.010. In the bottom panel, the data are shown relative to the standard 
model, shown as the solid line. The yellow diamonds are the result of binning the data (see table |B). 
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Fig. 16. — Our best-fit cosmological model in the S2 m -SiDB plane assuming a constant equation of state parameter w = —1. Our 
68%, 95%, and 99% c onfidence intervals are shown as solid blue (diagonal) contours. The constraints from baryon acoustic oscillations 
HEisenstein et al.ll2005h are shown as solid black (nearly- vertical) contours and the combined confidence intervals are shown as red contours. 
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Fig. 17. — Combining the CSP constraints with baryonic acoustic oscillations (Eiscnstcin ct al. 2005) and assuming Q.^ = 0. The CSP 
and BAO data combined are consistent with a value of w = — 1.05±0.13 (statistical) ±0.09 (systematic) and Q m = 0.27±0.02 (statistical). 
Our 68%, 9 5%, and 99% confiden ce intervals are shown as solid blue (banana-shaped) contours. The constraints from baryon acoustic 
oscillations (Eiscnstcin ct al. 2005]) are shown as solid black contours and the combined confidence intervals are shown as red contours. The 
1-D marginalized probabilities for each parameter are plotted as red lines on the axes. 
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Fig. 18. — Residuals in the i and B Hubble diagrams. The top panel shows reddening-free i magnitudes and the bottom panel shows 
the reddening-free B magnitudes. Both are computed using the (B — V) color. The rms dispersion is shown in the upper left. SN 04D2an 
is not shown as it has no V-band observation. If one includes SN 16442 (red square), the dispersion in the Hubble diagram increases to 
±0.15 at i, and ±0.18 at B as indicated in parentheses. The dashed line corrcponds to the best-fit cosmological model. 
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Fig. 19. — The sum (jj.) of the jerk (j) and curvature parameter (f2fc) as a function of the deceleration parameter (qo). The grey 
shading indicates the region where the luminosity distance expansion is valid, as described in the text. The best-fit values including both 
baryon acoustic oscillation and the CSP data are = 1.18 ± 0.44 (statistical) ± 0.27 (systematic) and q = —0.67 ± 0.13 (statistical) ± 
0.09 (systematic) at the 95% confidence level. The 1-D marginalized probabilities for each parameter are plotted as red lines on the axes. 
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Fig. 20. — The effect of a magnitude offset <5 m on the best-fit parameters, w, (Ide, Rvy B, Wo, and <tsn- The black circles are for 
i-band, while the red squares are for B-band. In almost all cases, the differences in the points are negligible. The lines are quadratic fits. 
The largest sensitivity to a magnitude offset is for w. For an offset of 0.025 magnitudes, the percentage change in w amounts to 7%. 
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Fig. 21. — The effect of a color offset 6 C on the best-fit parameters, as listed in Figure l20l The black circles are for i-band, while the 
red squares are for B-band. The lines are quadratic fits. The largest sensitivity to a color offset is again for w. For a color offset of 0.02 
magnitudes, the percentage change in w amounts to 5%. 
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Fig. 22. — The effect of extra variance in the color measurements cr c on the best-fit parameters, as listed in Figure l20l The black circles 
are for i-band, while the red squares are for B-band. The lines are quadratic fits. A variance in the color measurements impacts the 
determination of the reddening law or color term. A variance of ±0.03 magnitudes results in a percentage change in Ry of 10%. The 
impact on cosmology is very small. 
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FlG. 23. — The effect of an error in the color-term derived from local sequences of standards, as listed in Figure l20l The black circles arc 
for i-band, while the red squares are for B-band. The lines are quadratic fits. An error in the color term has the largest impact on w. For 
an error of 0.015, the percentage change in w amounts to 3%. 
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TABLE 1 
Filter nomenclature 



Name Eff. Wavelength 



Description 



Zero-point 



reference 



(Angstroms) 


Y 


10333.6 


Y-band on PANIC, Magellan I, LCO 


12 


6870 


J 


12440.3 


J-band on PANIC, Magellan I, LCO 


12 


8520 


B 


4905.0 


B-band on Swope, LCO 


13 


6672 


V 


4905.0 


V-band on Swope, LCO 


13 


7758 


i 


4905.0 


i-band on Swope, LCO 


13 


6873 


9m 


4905.0 


g-band on Megacam, CFHT 


15 


5363 


Trn. 


6282.6 


r-band on Megacam, CFHT 


14 


8061 


1 m 


7725.0 


j-band on Megacam, CFHT 


14 


5543 


Zm 


8901.6 


2-band on Megacam, CFHT 


13 


9996 


9a 


4718.9 


g-band on 2.5m, APO 


11 


2013 


r s 


6185.2 


r-band on 2.5m, APO 


14 


2157 


i s 


7499.7 


i-band on 2.5m, APO 


13 


7775 


Rim 


6660.2 


.R-band on 4m Blanco telescope, CTIO 


15 


1820 


I Am 


7973.8 


/-band on 4m Blanco telescope, CTIO 


14 


4580 



Persson et al. (_1998) 

55 

Hamuv et al. (2006') 



Astier et 



Smith et 



al. (2006 



Miknaitis et al. (2007) 
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TABLE 2 
CSP high-z SNe 



Name 


RA (2000) 
(d:m:s) 


DEC (2000) 
(li:m:s) 


Epochs 


z 


Cross-ID 


SNLS 03D4gl 


22:14:44.16 


-17:31:44.40 


3 


0.571 




SNLS 04Dloh 


02:25:02.38 


-04:14:10.68 


2 


0.59 




SNLS 04Dlpg 


02:27:04.15 


-04:10:31.33 


3 


0.515 




SNLS 04Dlrh 


02:27:47.16 


-04:15:13.68 


4 


0.435 




SNLS 04Dlsk 


02:24:22.56 


-04:21:13.32 


2 


0.6634 




SNLS 04D2an 


10:00:52.32 


02:02:28.68 


8 


0.62 




SNLS 05Dldn 


02:24:26.64 


-04:59:29.40 


1 


0.566 




SNLS 05Dlhk 


02:24:39.17 


-04:38:03.01 


2 


0.2631 




SNLS 05Dlhn 


02:24:36.26 


-04:10:54.95 


1 


0.1489 




SNLS 05Dlix 


02:24:19.94 


-04:40:11.75 


3 


0.49 




SNLS 05Dliy 


02:27:39.96 


-04:25:21.36 


1 


0.2478 




SNLS 05D2ah 


10:01:28.80 


01:51:46.08 


6 


0.184 




SNLS 05D2bt 


10:01:40.32 


02:33:57.96 


5 


0.679 




SNLS 05D2bv 


10:02:17.04 


02:14:26.16 


5 


0.474 




SNLS 05D2ck 


10:00:45.12 


02:34:22.08 


1 


0.698 




SNLS 05D2dw 


09:58:32.16 


02:01:56.28 


2 


0.417 




SNLS 05D2cb 


10:00:14.64 


02:24:26.64 


2 


0.5344 




SNLS 05D2mp 


09:59:08.64 


02:12:14.69 


3 


0.3537 




SNLS 05D4cw 


22:14:50.16 


-17:44:19.32 


3 


0.375 




SNLS 05D4fo 


22:15:20.88 


-17:16:05.16 


1 


0.373 




ESS dl49wcc4-ll 


02:10:53.98 


-04:25:49.80 


3 


0.342 


SN2003jy 


ESS el08wdd8-4 


02:30:09.00 


-09:04:35.76 


3 


0.469 


SN2003km 


SDSS 12855 


22:01:01.44 


00:42:58.47 


2 


0.165 




SDSS 13025 


22:46:16.08 


00:24:57.21 


3 


0.224 




SDSS 13835 


00:24:14.31 


-00:14:53.69 


2 


0.247 




SDSS 15287 


21:35:50.40 


-01:03:26.75 


1 


0.235 




SDSS 16442 


21:59:47.28 


-00:43:59.59 


3 


0.281 




SDSS 3241 


20:50:36.24 


-00:21:14.76 


5 


0.259 


SN2005gh 


SDSS 3331 


02:18:14.74 


00:47:47.76 


5 


0.208 


SN2005gc 


SDSS 4679 


01:26:06.79 


00:40:36.80 


1 


0.333 


SN2205gy 


SDSS 5183 


03:33:48.96 


00:42:33.70 


6 


0.384 


SN2005gq 


SDSS 5549 


00:13:00.13 


00:14:53.70 


5 


0.119 


SN2005hx 


SDSS 6699 


21:31:15.60 


-01:03:25.16 


1 


0.311 


SN2005ik 


SDSS 7243 


21:52:18.96 


00:28:19.09 


3 


0.204 


SN2005jm 


SDSS 7512 


03:28:21.67 


-00:19:34.10 


2 


0.22 


SN2005jo 
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TABLE 3 
CSP high-z Photometry 



Name 


MJD a 


Y c 




K Y ,i b 


Jc 


(J J 


K J,i c 




fdavs) 


( rnatr) 


fmae'l 


( mac) 


(mag) 


(mae) 


(mae*) 


SNLS 03D4gl 


52947.0 


22.57 


0.09 


-0.97 


22 


84 


0.15 


-1.37 




52947.0 


22.54 


0.05 


-0.91 


23 


03 


0.13 


-1.40 




52952.0 


22.64 


0.05 


-0.90 


22 


80 


0.08 


-1.38 


SNLS 04Dloh 


53302.1 








22 


81 


0.09 


-1.39 




53314.2 








22 


78 


0.10 


-1.37 


SNLS 04Dlpg 


53328.2 








22 


72 


0.08 


-1.46 




53338.1 








22 


95 


0.08 


-1.48 




53344.2 








23 


28 


0.19 


-1.59 


SNLS 04Dlrh 


53356.1 


22.44 


0.06 


-1.21 


22 


26 


0.07 


-1.34 




53356.1 


22.60 


0.09 


-1.23 


22 


53 


0.11 


-1.46 


SNLS 04Dlsk 


53356.1 








23 


33 


0.19 


-1.50 




53363.0 








23 


78 


0.33 


-1.51 


SNLS 04D2an 


53027.2 


22.75 


0.05 


-1.11 


22 


92 


0.07 


-1.41 




53027.2 


22.63 


0.05 


-1.08 


23 


23 


0.11 


-1.40 




53034.1 


22.84 


0.10 


-1.03 


23 


07 


0.14 


-1.39 




53034.1 


22.67 


0.05 


-0.95 


23 


58 


0.17 


-1.39 




53040.2 


23.11 


0.07 


-0.91 


23 


64 


0.17 


-1.41 




53040.2 


23.27 


0.12 


-0.90 


23 


52 


0.22 


-1.41 




53045.2 


23.32 


0.07 


-1.03 


23 


38 


0.14 


-1.41 




53045.2 


23.33 


0.10 


-1.06 


23 


73 


0.20 


-1.41 


SNLS 05Dldn 


53641.3 








22 


76 


0.13 


-1.35 




53654.3 








22 


64 


0.10 


-1.35 




53659.2 








22 


94 


0.13 


-1.36 




53666.2 








23 


11 


0.19 


-1.41 


SNLS 05Dlhk 


53720.0 
53724.1 


21.10 
21.13 


0.03 
0.03 


-1.02 
-0.98 










SNLS 05Dlhn 


53721.0 


20.63 


0.03 


-0.97 










SNLS 05Dlix 


53720.1 








22 


07 


0.09 


-1.54 




53722.1 








22 


27 


0.09 


-1.54 




53724.0 








22 


04 


0.09 


-1.51 


SNLS 05Dliy 


53721.1 


21.65 


0.02 


-0.92 










SNLS 05D2ah 


53390.2 
53390.2 


20.84 
20.93 


0.03 
0.02 


-0.86 
-0.87 


21 


46 


0.05 


-0.39 




53391.2 








21 


71 


0.06 


-0.38 




53392.2 


20.93 


0.02 


-0.92 


21 


71 


0.08 


-0.37 




53393.3 


21.02 


0.03 


-0.95 


21 


77 


0.07 


-0.37 




53393.3 


20.93 


0.02 


-1.16 


21 


66 


0.04 


-0.37 


SNLS 05D2bt 


53405.1 








22 


98 


0.09 


-1.54 




53413.2 








23 


28 


0.24 


-1.55 




53418.1 








23 


07 


0.34 


-1.56 




53424.1 








23 


47 


0.28 


-1.54 




53430.1 








23 


09 


0.20 


-1.54 


SNLS 05D2bv 


53405.3 








22 


21 


0.04 


-1.35 




53413.2 








22 


64 


0.13 


-1.36 




53418.2 








22 


58 


0.08 


-1.41 




53424.2 








22 


76 


0.15 


-1.58 




53430.1 








22 


90 


0.14 


-1.60 


SNLS 05D2ck 


53418.3 








23 


28 


0.25 


-1.54 




53424.3 








23 


29 


0.30 


-1.58 




53430.2 








23 


51 


0.23 


-1.59 




53444.1 








23 


61 


0.20 


-1.55 


SNLS 05D2dw 


53447.1 


21.99 


0.07 


-1.16 


22 


20 


0.09 


-1.31 




53447.1 


22.20 


0.06 


-1.14 


22 


00 


0.06 


-1.30 


SNLS 05D2eb 


53447.1 
53447.1 


22.23 
22.49 


0.12 
0.04 


-0.82 
-0.78 


22 


10 


0.15 


-1.51 


SNLS 05D2mp 


53720.3 
53722.2 
53724.3 


22.44 
22.55 
22.54 


0.04 
0.04 
0.04 


-0.94 
-0.95 
-0.96 










SNLS 05D4cw 


53591.1 
53596.1 
53608.1 


22.35 
22.66 
22.55 


0.07 
0.08 
0.08 


-1.01 
-1.03 
-1.05 










SNLS 05D4fo 


53648.1 


22.14 


0.45 


-1.00 


22 


07 


0.15 


-1.35 




53648.1 


22.32 


0.07 


-1.00 


21 


64 


0.06 


-1.50 




53654.1 


22.50 


0.12 


-1.00 


21 


79 


0.11 


-1.48 




53654.1 


22.76 


0.10 


-1.00 


21 


96 


0.11 


-1.56 


ESS dl49wcc4-ll 


52947.3 
52952.2 
52956.2 


22.01 
21.92 
21.97 


0.04 
0.03 
0.03 


-0.91 
-0.92 
-0.91 










ESS el08wdd8-4 


52980.1 


22.42 


0.04 


-1.29 


22 


25 


0.08 


-1.38 




52980.1 


22.24 


0.03 


-1.30 


22 


28 


0.09 


-1.39 




52982.1 


22.20 


0.03 


-1.31 


22 


28 


0.08 


-1.39 


SDSS 12855 


54000.1 
54008.0 


20.61 
20.92 


0.02 
0.02 


-0.89 
-1.10 










SDSS 13025 


53996.1 


21.02 


0.03 


-1.00 
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TABLE 3 — Continued 



Name 


MJD a 

(davs) 


Y c 
( mae) 


<7y 
fmae:^ 


K Y r h Jc CFJ KjT C 
(masl (mag) (maer) fmae;') 




54000.1 


21.08 


0.04 


-0.99 




54008.1 


21.58 


0.04 


-1.20 


SDSS 13835 


54008.1 


21.14 


0.01 


-0.90 




54021.2 


21.40 


0.02 


-1.01 


SDSS 15287 


54029.0 


21.13 


0.02 


-0.94 


SDSS 16442 


54060.0 


21.30 


0.04 


-0.86 




54065.0 


21.46 


0.03 


-0.80 




54070.0 


21.79 


0.04 


-0.80 


SDSS 3241 


53641.0 


21.31 


0.02 


-0.86 




53648.0 


21.27 


0.02 


-0.93 




53654.0 


21.54 


0.02 


-0.90 




53659.0 


21.79 


0.02 


-0.95 




53666.0 


21.98 


0.03 


-1.05 


SDSS 3331 


53641.2 


20.80 


0.03 


-0.88 




53648.3 


20.50 


0.04 


-1.06 




53655.2 


20.77 


0.03 


-1.03 




53659.2 


20.82 


0.03 


-1.10 




53676.2 


20.89 


0.04 


-1.32 


SDSS 4679 


53655.1 


21.99 


0.02 


-0.89 




53659.1 


22.02 


0.03 


-0.88 




53666.2 


22.35 


0.04 


-0.89 




53676.2 


22.37 


0.03 


-0.91 


SDSS 5183 


53654.3 


22.39 


0.06 


-1.03 




53655.3 


22.46 


0.05 


-1.04 




53666.2 


23.22 


0.09 


-1.04 




53676.3 


22.94 


0.07 


-1.06 




53720.2 


23.68 


0.14 


-1.09 




53721.2 


23.72 


0.16 


-1.09 


SDSS 5549 


53666.1 


19.88 


0.04 


-0.56 




53676.1 


20.20 


0.04 


-0.65 




53682.1 


20.03 


0.11 


-0.93 




53684.1 


20.18 


0.03 


-0.95 




53689.1 


20.11 


0.04 


-0.97 


SDSS 6699 


53676.0 


21.72 


0.02 


-0.90 




53682.0 


21.87 


0.05 


-0.85 




53684.0 


22.16 


0.04 


-0.84 




53689.0 


22.42 


0.05 


-0.85 


SDSS 7243 


53682.1 


20.82 


0.02 


-0.83 




53684.1 


20.87 


0.02 


-0.88 




53689.0 


20.94 


0.01 


-0.88 


SDSS 7512 


53682.3 


21.21 


0.10 


-0.89 




53684.3 


21.25 


0.07 


-0.89 



a MJD = JD - 2400000.5. 

k K-correction from observed filter Y c to rest-frame filter I. 

K-correction from observed filter J c to rest-frame filter I, except for SNLS05D2ah, for which the K-correction is to rest-frame Y. 
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TABLE 4 

Derived Light-curve Parameters for low-z SNe 



Name 




z 




a (oY 


J max \o ) 






°) 










(mag 


) 


(days) 




(ma 


I) 


SN2004ef 





03097 


35 


A Qo/n 


ioy ) 


264.8(0.0) 


1 


389 





007) 


SN2004eo 





01569 


33 


oy ( i^u 


ioy ) 


278.9(0.1) 


1 


366 





012) 


SN2004ey 


o 


01578 


3 1 


no7/ n 
UZ ( \\) 


iyu ) 


304.6(0.0) 





954 





006) 


fiN2004ffs 


o 


02663 


35 


A Q.A ( C\ 

4o4|^U 


1 A Q~\ 

14z ) 


356 3?D f)^ 


1 


550 





006) 


SN2004gu 


o 


04583 


36 


A QA f n 
4o4|^U 


izy ) 


362.1(0.2) 


o 


758 





011) 


SN2005ag 


o 


07937 


37 


d A Q (c\ 
u4o(U 


1Z1 J 


414.1(0.1) 


o 


889 





007) 


SN2005al 


o 


01239 


34 


uyz\u 


ZUZ ) 


430 6f0 1 ) 


1 


243 





010) 


SN2005ol 


o 


01490 


33 


yo ( \\) 


toy j 


647.0(0.1) 


1 


299 





014) 


SN2005en 


o 


02896 


35 


489(0 


141) 


654 5^0 I) 


o 


778 





008) 


SN2005hc 


o 


04591 


36 


582(0 


128) 


667 4f0 1 1 


o 


844 





008) 


SN2005hj 





05797 


36 


976(0 


127) 


674.0(0.2) 





739 





017) 


SN2005iq 





03402 


35 


891(0 


136) 


687.9(0.1) 


1 


230 





016) 


SN2005ir 





07631 


37 


589(0 


127) 


685.0(0.2) 





875 





025) 


SN2005kc 





01511 


33 


871(0 


197) 


698.3(0.0) 


1 


150 





016) 


SN2005ki 





01919 


34 


615(0 


159) 


705.9(0.0) 


1 


381 





007) 


SN2005M 





02200 


35 


052(0 


151) 


405.9(0.0) 





799 





003) 


SN2005na 





02630 


35 


118(0 


145) 


740.3(0.2) 


1 


005 





013) 


SN2006ax 





01673 


34 


324(0 


169) 


827.5(0.1) 





949 





008) 


SN2006bh 





01084 


33 


329(0 


238) 


833.6(0.0) 


1 


387 





008) 


SN2006gt 





04474 


36 


441(0 


130) 


1003.0(0.1) 


1 


675 





012) 


SN2006py 





05786 


36 


839(0 


137) 


1071.0(0.3) 


1 


016 





049) 



(imo) (<7) d (Bmax-M W E{B - V) host (<r)° A TO T C 





(mag 


) 


(mag) 


(mag) 


(mag) 


17 


273 





010) 


0.123 


(0.007) 


0.157 


(0.018) 





214 


15 


435 





016) 


0.095 


(0.015) 


0.129 


(0.021) 





249 


15 


459 





011) 


-0.070 


(0.008) 


-0.003 


(0.015) 





163 


1 7 


391 





uuy / 


0.202 


(0.009) 


0.223 


(0.022) 





247 


18 


052 





020) 


0.145 


(0.021) 


0.253 


(0.027) 





268 


19 


123 





011) 


0.004 


(0.007) 


0.085 


(0.016) 





128 




609 





m ^ 

U-LO f 


-0.089 


(0.009) 


-0.057 


(0.016) 





012 


15 


538 





018) 


-0.055 


(0.014) 


-0.025 


(0.020) 





112 


ID 


962 





U±0 ) 


0.039 


(0.010) 


0.135 


(0.020) 





215 


18 


050 





014) 


0.009 


(0.008) 


0.095 


(0.017) 





129 


18 


462 





026) 


0.067 


(0.014) 


0.170 


(0.023) 





207 


17 


457 





019) 


-0.033 


(0.011) 


0.005 


(0.017) 





032 


19 


086 





037) 


0.027 


(0.016) 


0.112 


(0.022) 





141 


15 


632 





020) 


0.200 


(0.014) 


0.267 


(0.019) 





409 


16 


237 





010) 


-0.031 


(0.010) 


-0.010 


(0.019) 





029 


16 


518 





006) 


0.024 


(0.004) 


0.117 


(0.017) 





146 


16 


618 





020) 


-0.018 


(0.015) 


0.048 


(0.019) 





137 


15 


763 





014) 


-0.060 


(0.011) 


0.009 


(0.017) 





068 


14 


977 





010) 


-0.007 


(0.009) 


0.016 


(0.018) 





047 


18 


415 





020) 


0.225 


(0.020) 


0.234 


(0.031) 





264 


18 


421 





060) 


0.070 


(0.021) 


0.142 


(0.025) 





205 



a Distance modulus (h = 0.72). 

k Time of maximum for rest-frame B light-curve (JD - 2400000.5). 
c Decline rate parameter. 

^ Galactic reddening from lSchleeel et alj l| 19981 ). 

Host galaxy reddening assuming IPhillips etaLl ( fT999f) colors. 

r 

Total (galactic + host galaxy) absorption in NIR toward SN assuming Ry = 1.7. 
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TABLE 5 
Derived Light-curve Parameters 



Name 




z 




(mag) 


Tmax \& ) 

(days) 


1) 




(ma 


r) 


H 


max ) 

(mag 


<7) d 
) 


(Bmax — Vmax) (o") d 

(mag) 


E(B - 


V)h„st (cr) c 
(mag) 


(mag) 


SNLS 


03D4gl 





571 


42 


765 





170) 


52954 


1 





6) 





719 





151) 


24 


192 





120) 


0.020 (0.029) 


0.122 


;0.038) 


0.136 


SNLS 


04Dloh 





59 


42 


680 





180) 


53306 


9 





4) 





946 





069) 


21 


070 





087) 


-0.026 (0.056) 


0.046 


;0.058) 


0.067 


SNLS 


04Dlpg 





515 


42 


353 





148) 


53325 


5 
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a Distance modulus (h = 0.72). 

k Time of maximum for rest-frame B light-curve (JD - 2400000.5). 
c Decline rate parameter. 

^ Galactic reddening from [Schlcgcl ct al.l 1 119981 ). 

Host galaxy reddening assuming IPhillips etaTI JT999D colors. 

r 

Total (galactic + host galaxy) absorption in NIR toward SN assuming Ry = 1.7. 
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TABLE 6 
Binned Hubble Diagram 



< z > 


N 


< n> 


&( u) 

Vr*7 


0.01 


6 


33.906 


0.018 


0.02 


3 


34.635 


0.049 


0.03 


5 


35.505 


0.028 


0.05 


5 


36.619 


0.032 


0.08 


2 


37.595 


0.045 


0.13 


2 


38.930 


0.117 


0.17 


2 


39.627 


0.099 


0.24 


10 


40.357 


0.054 


0.35 


7 


41.397 


0.057 


0.46 


5 


41.984 


0.065 


0.56 


5 


42.543 


0.075 


0.68 


3 


43.023 


0.107 
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TABLE 7 

COSMOLOGICAL RESULTS 



R v Band Q m w rms 

1.74 Wty 0.27 ± 0.02 -1.05 ±0.13 0.13 

Wg v 0.27 ±0.02 -1.08 ±0.14 0.15 

3.1 Wty 0.26 ±0.02 -1.20 ±0.13 0.20 

WS V 0.25 ±0.03 -1.24 ±0.16 0.24 
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TABLE 8 
Systematic Errors 









= 


UI 


-1 




j =constant 








nuisance 






Systematic 


Error 


dw 
da 


Aw 


da 


An DE 


dqo 

da 


Ago 


•in, 

da 




da 


AR V 


rib 
da 


A6 


da 


AW 


Mag. offset (8m) 


0.025 


-2.69 


0.067 


2.73 


0.068 


-2.75 


0.069 


8.10 


0.203 


0.04 


0.001 


0.002 


0.000 


-0.25 


0.006 


Color offset (8c) 


0.02 


2.53 


0.051 


-2.58 


0.052 


2.68 


0.054 


-7.77 


0.155 


0.38 


0.008 


-0.03 


0.001 


0.21 


0.004 


Color Term (8(CT)) 


0.015 


1.8 


0.027 


-1.73 


0.026 


1.83 


0.027 


-5.82 


0.087 


-0.67 


0.010 


-0.22 


0.003 


0.71 


0.011 


Amis variance (cta) 


0.05 


-0.007 


0.000 


-0.01 


0.001 


-0.04 


0.002 


0.14 


0.007 


0.09 


0.005 


-0.06 


0.003 


-0.03 


0.002 


Color variance (a c ) 


0.03 


0.03 


0.001 


0.52 


0.016 


-0.173 


0.005 


0.53 


0.016 




0.100 




0.011 




0.004 


Total 






0.088 




0.091 




0.092 




0.270 




0.097 




0.012 




0.014 
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TABLE 9 

Summary of Cosmological Parameters 



Constraints Cosmological Parameters 

w = -1 Q m = 0.27 ± 0.02 (statistical) 

Combine with BAO n DE = 0.76 ± 0.13 (statistical) ± 0.09 (systematic) 

n k = n m = 0.27 ± 0.03 (statistical) 

Combine with BAO w = —1.05 ± 0.13 (statistical) ±0.09 (systematic) 



Constant j q = —0.67 ±0.13 (statistical) ±0.09 (systematic) 

Combine with BAO j k = 1.18 ± 0.44 (statistical) ± 0.27 (systematic) 



